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Abstract 

In this paper we prove that the two dimensional superintegrable 
systems with quadratic integrals of motion on a manifold can be clas- 
sified by using the Poisson algebra of the integrals of motion. There 
are six general fundamental classes of superintegrable systems. Ana- 
lytic formulas for the involved integrals are calculated in all the cases. 
All the known superintegrable systems are classified as special cases 
of these six general classes. 



Running title: Classification of two dimensional superintegrable systems 
PACS Numbers: 03.65.Fd; 02.10.Tq; 45.20. Jj; 



2 



I Introduction 



In classical mechanics, an integrable system on a manifold of N dimensions, 
is a system which has N integrals of motion in evolution. Superintegrable (or 
maximally integrable) system is a system possessing the maximum number 
of constants of motion, i.e. 2N — 1 integrals of motion. 

The simplest case is the two dimensional superintegrable problems with 
integrals of motion, which are linear and quadratic functions of momenta. 
The investigation of such superintegrable systems on a two dimensional man- 
ifold is a quite old problem, dated on 19th century. Initially the problem was 
formulated as a geometry problem. The challenge was to find two dimen- 
sional manifolds whose the geodesies are curves which possess additional 
integrals than the free Hamiltonian. This problem was studied in the four 
volume treatise of Darboux Legons sur la Theorie Generale des Surfaces pQ. 
The main result of this study is that, there are five classes of general forms 
of metrics, whose the geodesies have three integrals of motion (the Hamilto- 
nian and two additional functionally independent integrals). These metrics 
are called "formes essentielles" and they depend on four parameters. All the 
metrics having more than two integrals of motion can be obtained as partial 
cases of these "formes essentielles" by choosing appropriate values of the four 
parameters. The five classes of metrics are tabulated in "Tableau VII" by 
Koenigs[U vol IV p. 385] 

In Classical Mechanics language, Darboux and Koenigs results can be 
translated as searching manifolds, where the free Hamiltonian accepts quadratic 
integrals of motion. The evolution of this problem is to find superintegrable 
systems, whose Hamiltonian is the free Hamiltonian plus a potential and 
these systems possess additional quadratic integrals of motion. 

The simplest integrable and superintegrable systems are the systems de- 
fined on the real plane. A comprehensive review of the real two-dimensional 
integrable classical systems on the plane is given in Ref [21 [3]. The com- 
plex superintegrable systems with quadratic moments on a flat space were 
recently catalogued by Kalnins, Miller, Pogosyan et al [U El El E] • The flat 
space is a two dimensional manifold with curvature zero. The Drach poten- 
tials are also systems defined on a manifold with curvature zero. The Drach 
systems with quadratic integrals of motion were investigated by Ranada [8] . 
The superintegrable systems on the hyperbolic plane were studied in [91 [10] . 
These systems were studied separately, while they are connected by obvi- 
ous coordinate transformations. Therefore they are naturally connected, i.e. 
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there is a common classification scheme of all these systems [HI [TQ] . 

The case of non flat space is under current intensive investigation. The 
superintegrable systems on the sphere were studied in |TT] and they 
were classified in [6]. The sphere is a special case of manifold with constant 
curvature. In refs (9J [10] the superintegrable systems on the sphere and on 
the hyperbolic plane were studied using an unified formulation. 

In the case of manifolds of non-constant curvature, the known examples 
of superintegrable systems are those which are defined on manifolds which 
are surfaces of revolution. The corresponding problem of the geodesies with 
three quadratic integrals of motion on surfaces of revolution was treated 
by Koenigs in [Tj n° 5, vol IV, p. 377]. Recently Kalnins and collaborators 
classified the superintegrable systems on a surface of revolution [TJl [32] using 
the manifolds which were provided by Koenigs. 

In two recent papers Kalnins, Kress and Miller jHJ [15] give a compre- 
hensive study of the two-dimensional superintegrable sytems. In [15], they 
prove that the general Koenigs essential forms of metrics correspond to the 
most general forms of superintegrable systems. Also they have shown that 
every two-dimensional superintegrable system is Stackel equivalent to a two- 
dimensional non degenerate superintegrable system on a constant curvature 
space. 

Kress [16] in collaboration with Kalnins and Miller studied the Stackel 
equivalence classes of the superintegrable systems on the spaces of constant 
curvature, and they have shown that there are six equivalence classes. The 
general classes of Koenigs classification given by the Table VII in [lj are 
five. That means that there is a sixth class which should be added in the 
Koenigs classification scheme. In this paper we investigate this sixth class, 
which completes the Koenigs classification. This class is the nondegenerate 
superintegrable system generated by the case VI 6 of Koenigs. 

An interesting question is, whether there could be a general classification 
scheme of superintegrable systems with quadratic integrals of motion, which 
contains all the equivalence classes of superintegrable systems on a manifold 
with constant curvature and the general classes of manifolds , which were 
introduced by Koenigs. The classification schemes are based on the Darboux 
relations derived for the invariants, which are defined on a specific manifold. 
In this paper we propose a classification scheme based on the properties of 
the Poisson algebra of the integrals of motion. Then we show that there is in- 
deed such a classification scheme, which determines the supporting manifold 
metric. We must notice that the Kress [T6] equivalence classes are derived by 
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classifying the Poisson algebra of integrals of motion. The proposed classes 
of superintegrable systems in this paper correspond to the equivalent classes 
studied in ref [16] . Analytic formulas for the metrics of the permitted mani- 
folds, the potentials and the integrals of motion are calculated. 

This paper is organized as follows: In section [III the general form of the 
integrable two dimensional system with one quadratic integral of motion is 
derived. The results of this section correspond to the Darboux treatise para- 
graphs [H n° 593, vol. Ill, p. 30], [TJ n° 593, vol. Ill, p. 31], but we give a brief 
modern derivation of the formulas including the potentials in our discussion. 
These formulas will be used in the following sections. The carrying manifold 
is a Liouville or a Lie surface. Therefore there are two classes of integrable 
systems. In a specific coordinate system, which is called Liouville (or Lie) 
system, the analytic expressions of the potential and the integrals of motion 
are given and the action is calculated. In Section II III the Poisson algebra of 
the integrals of a superintegrable two dimensional system is discussed. This 
algebra is a quadratic algebra, the coefficients of the quadratic terms are 
characteristic of the carrying manifold. In Section HVl we prove that the coef- 
ficients of the Poisson algebra impose the classification of the superintegrable 
systems with two quadratic integrals in six fundamental classes. The method 
of analytic calculation of form of the permitted carrying manifolds, the poten- 
tials and the integrals of motion are discussed. In this section we prove that 
the general form of the superintegrable potential can be written as a fraction 
V = w(x, y)/g(x, y) and the two functions w(x, y) and the metric g(x, y) are 
two solutions of the same partial differential equation. The existence of the 
Poisson algebra was assumed as obvious by several authors [I] [I], [T7]-[29j. 
In Appendix [A] we give a proof of the existence of the Poisson algebra, for 
two dimensional superintegrable systems with quadratic integrals of motion. 
In Section [V] the analytic formulas of the manifolds and integrals are given 
for all the six fundamental classes of superintegrable systems. From these 
analytic formulas we can show that there are new superintegrable systems, 
because they are defined on manifolds which have not constant curvature and 
are not surfaces of revolution. In Section [VI] the superintegrable systems cor- 
responding to the Koenigs essential forms of Table VII are given. In Section 
IVHI the superintegrable systems on a surface of revolution are studied. We 
find that there is a new system which was not revealed by the other classi- 
fication schemes. In Section IVIIII the superintegrable systems on a manifold 
with curvature zero are studied and in Section ITXl the systems on a manifold 
with constant curvature are listed. In Section [X] the systems with a linear 
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and a quadratic integral of motion are discussed. Finally, in Section IXII the 
results of the paper are summarized. 



II Integrable systems on a two dimensional 
manifold 

Let us consider an integrable system defined on a two dimensional manifold 
with metric: 

ds 2 = E(u, v)du 2 + 2F(u, v)dudv + G(u, v)dv 2 

There is a conformal coordinate system where the metric can be written as: 

ds 2 = g(x,y)dxdy (1) 

The passage from the original coordinate system (u, v) to the conformal one 
(x, y) can be realized by using the Beltrami partial differential equation. We 
must notice that the choice of the conformal coordinate system is not unique, 
i.e. there are several conformal coordinate systems for a given metric, these 
systems are conformally equivalent. 

In a conformal coordinate system the general form of the Hamiltonian is: 

where the Hamiltonian is a quadratic form of the momenta. 

Let us consider an integral of motion which is quadratic in momenta. The 
most general form can be written as: 

/ = A(x, y) V 2 x + B(x, y) V 2 y - 2p xPy ^\ + Q(x, y) (3) 

By definition the Poisson bracket between the Hamiltonian and the integral 
of motion is zero: 

_ dl_dH_ _ dl_dH_ + d/aff _ dl_dH_ _ Q 

p dx dp x dp x dx dy dp y dp y dy 

The above equality implies restrictions on the functions involved in equations 
(j2J) and (J3j). The left hand side of the above equation is an odd function of 
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cubic order in the momenta. The coefficients of p^ and p\ must be zero: 

= => A = A(x) 

dB 



dA 
dy 



(5) 



dx * B = B(y) 
The coefficients of and in (TjJ must be indeed zero: 

= A{x) dx- + 2 A {X) = a{x) dx (6) 

Tx = B{v) Wy + 2 B iy) = KV) dy (7) 

where 

A(x) = a 2 (x), B(y) = b 2 (y) 

The partial x-derivative of the right hand side in ([2D is equal to the y- 
derivative of the right hand side in equation (17j) , therefore: 

(A"(x) - B"(y))g(x, y) + 3A'(x)?jL - 3B'(y)?f + 2A(x)^ - 2S(y) = 

ox ay dx z dy z 

(8) 

or 

did \ d ( d \ 

dx~ [^di WWW' = ~dy~ [ b ^dy~ V ^ J (9) 

The coefficients of p x and p y in must be zero: 



dy dx dy 



(10) 



The above relations imply: 



g(x,y) {2A(x)^ - 2B{y)% + 3A'{x)% - W{y)%) + 
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At this point we must to distinguish two cases. In the first case A(x) and 
B(y) are both different from zero, whereas in the second case B(y) is assumed 
to be zero. 

Class I: A(x)B(y) f 

Following the method given in Koenigs original paper, we can choose a new 
coordinate system 

dx f dy 



£ = f ; and ri 



where the associated momenta are: 

PS, = \j A i x )Px, and p v = ^B(y)p y 

In this case the metric is written 

ds 2 = g(£, rj) d£ dr], where g(£, rf) = g(x, y) sj A{x) y ' B(y) 

In these new coordinate system (£, rj) all the above equations are considerably 
simplified. We can easily show that the formulas are the same by replacing 
x — > £, y — > rj and fixing A(x) = 1 and B(y) = 1. For simplicity reasons 
we omit the hat on the metric g(£,,T]) — > g(C,v)- Equations (|2J) and are 
written: 

2 o P(Z,V) , „ft X (12) 



We call these specific coordinates (£, r/) Liouville coordinates . In Liouville 
coordinates the Hamiltonian H and the integral / are written as it was given 
in the above equation (|T2l) . 

In the Liouville coordinates, equation (jUj) is considerably simplified: 

d 2 g d 2 g _ 

the general solution is: 

9&ri) = F{Z + v) + G{Z-7i) 

where F(u) and G(v) are arbitrary functions. The above metric character- 
izes a surface which is called Liouville surface in the geometry textbooks. 
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Therefore we have shown that a system is integrable in Class I only when 
the surface is a Liouville surface. The surfaces of constant curvature or the 
rotation surfaces are Liouville surfaces, but there are Liouville surfaces which 
have not constant curvature and they are not rotation surfaces. 
The function /?(£,7?) can be calculated from equations (jUJ) 

d/3 dg d(3 dg 
<9£ dr]'' drj <9£ 

then 

P(Z, V ) = F(£ + t 1 )-G(S-t)) (13) 
The potential V(£, 77) in Liouville coordinates is the solution of equation ( TTTT) 



(F« + ,) + C« - „)) - 0) + 
+ 2F<« + ,)(f -f)+2G'« + ,)(f + f)=0 

The general solution of the above equation in Liouville coordinates is 

The functions /(m) and g(t>) are arbitrary functions. The function Q(£,rj) is 
determined from equations ffTUl) and the solution is easily calculated: 

Q({ '" ) = 4 F(£ + ,)+G(£-,) (15) 

Usually is more convenient the use of the coordinates u, t>, which are defined 
by: 

t. . ■ Pu - iPv , P« + ip« 
^ = M + zu, 7] = u-iv, p ( = , and p v = 

The Hamiltonian H and the integral can be written as: 

„_ Pu 2 +Pv 2 + 4:(f(u) + g(v)) 
4(F(u) + G(v)) 

T Pu 2 G(v) -p v 2 F(u) , A f(u)G{v) - g(v)F(u) 



+ 4 



F(u) + G(v) F(u) + G(v) 
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The above formula has been investigated in a different context in [30] . In 
these coordinates the action S(u, v) satisfy the following equations: 

E = H(u, v,p u ,p v ), J — I(u,v,p u ,p v ) , . 

„ _ OS _ dS l iC) J 

Pu du , Pv dv 

and we can find the action S(u, v) by separation of variables. 

S = -Et + J \JAE F{u) + J-4f(u) du + J ^4EG(v) -J-4g(v) dv 

CLass II: B(y) = 

We can choose a new coordinate system 

_ r dx _ 

the associated momenta are: 

PZ = \jA(x)p x , and p v = p y 

In this case the metric is written 

ds 2 = g(£, rj) d£ drj, where g(£, 77) = g(x, y) \j A(x) 

We can easily show that the formulas are the same by replacing x — > £, y — > r] 
and fixing A(x) = 1 and B(y) = 0. For simplicity reasons we omit the hat 
on the metric g{C,v) ~^ di.^ 7 ])- Equations ^ and ([3]) are written: 

, ^ \ (17) 



I = p\- 2p i p 1) — + Q(f , 77) 



We call these specific coordinates (£, 77) Lie coordinates . In Lie coordinates 
the Hamiltonian H and the integral / are written as it is given in the above 
equation ( jTTl) . 

In Lie coordinates equation flU]) is written: 

^=0 
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The general solution is: 



g(Z,v) = F( V )Z + G(ri) 

where F(rj) and G(rj) are arbitrary functions. The above metric characterizes 
a surface which will be called Lie surface. Therefore we have shown that a 
system is integrable in Case 2 only when the surface is a Lie surface. 
Equations imply: 



df3 _dg 1 
drj <9£ 



dp 







g(Z,v)=F(ri)Z + G(7i) 
P&V) = f F(ri) dr, 



(18) 



Equation lull) is written: 



and the general solution of the above equation is: 



F(r))£ + G(r)) 



(19) 



The functions F(r]), G(r)), f{rj) and g{rj) are arbitrary functions. In this 
case the solution of the system of equations (ITU]) is given by: 



,(f(v)£ + g(v))f Hv) df] 



F(v) £ + G{rj) 

The action integral in this case can be easily calculated 



+ 2 J f(v) dr, 



(20) 



S = ^J-2[ I f( V )d V -Ej F( V )d V 
where 



g{r,) - EG(rj) 



J-2(f f(r))d V - Ef F(r))dr)) 



The above findings are summarized in the following Proposition. 
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Proposition 1 A Hamiltonian H quadratic in momenta, which is defined 
on a two dimensional manifold possesses an integral of motion I quadratic 
in momenta, only in two cases: 
Class I: 

The manifold is a Liouville surface, i.e. there is a coordinate system (£, rf) 
where the metric can be written: 

ds 2 =g(^r])d^ and g(£, 77) = F(£ + rj) + G(f - 77) 

H = P±Pv /(£ + ?7) + - ?7) 



+ V) + G(£ - v) F(£ + 77) + G(£ - r?) 
and simultaneously the integral of motion is written 

Class II: 



The manifold is a Lie surface, i.e. there is a coordinate system (£, 77) where 
the metric can be written: 

ds 2 = g(£,rj)d£dr) and g(£,rj) = F(rj) £ + G(rj) 

H _ PzPy f(v)€ + g(v) 



F( V )^ + G( V ) F( V )C + G(rj) 
and simultaneously 
, 2 , ! F{r,)d v (/(,)) £ + g(r,)) { F(r,) dr, ,„/■-,. 

The above specific choice of coordinate system (£, 77) will be called Lie coor- 
dinate system. 

The integrable systems which belong in Class I and II are well known inte- 
grable systems see Ref. [12], the supporting manifold is a Liouville or a Lie 
surface. 

We must notice that the Hamiltonian and the integral of an integrable 
system determine uniquely the Liouville (or Lie) coordinate system. There- 
fore the use of this privileged system is imposed by the notion of integrability. 
In the next sections we shall work in this special coordinate system, which is 
denoted exclusively by the coordinates £ and 77. In Class I integrable systems 
the system is separable, while in Class II systems there is no separation of 
variables in general. 
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Ill Poisson algebra of superintegrable systems 
with two quadratic integrals of motion 

If a system is superintegrable on a two dimensional manifold, that means 
that there are three functionally independent integrals of motion H, A and 
B. In this section, we assume that these integrals of motion are quadratic 
functions of the momenta and there are no other integrals of motion, which 
are linear functions of the momenta. Regarding the Hamiltonian H and the 
first integral A, we can choose the Liouville coordinate system and in this 
system: 

As we have shown in the previous section the system is integrable with a 
square integral of motion in two cases. The integral of motion A in the 
Liouville coordinate system is written: 

a = P l + bpl- 2 P ,p^|^ + e(£, ?7 ) 

where 

1 in Class I, where rj) = F(£ + rj) + G(£ — rj) (Liouville system) 

< 

in Class II, where g(£,r)) = F(rj) £ + G(r)) (Lie system) 



b= I 



The integral of motion B is assumed to be indeed a quadratic function of 
momenta, thus the general form in Liouville coordinates is 

B = A(0?l + B(rj)pl - 2p^f|g + Q(£ rj) 

By definition the following relations are satisfied: 

{H,A} p = {H,B} p = (21) 

From the integrals of motion A and B, we can construct the integral of 
motion: 

C = {A,B} P (22) 

The integral of motion C is not a new independent integral of motion, which 
is a cubic function of the momenta. The integral C is not functionally inde- 
pendent from the integrals H, A and B as it will be shown later. The fact 
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that, the integral C is a cubic function of momenta, implies the impossibility 
of expressing C as a polynomial function of the other integrals, which are 
quadratic functions of momenta. We shall prove that the square of this cu- 
bic polynomial is indeed a cubic combination of the integrals. Starting from 
the integral of motion C, we can construct the (non independent) integrals 
{A, C} p and {B, C} p . These integrals are quartic functions of the momenta, 
i.e. functions of fourth order. In Appendix |A] we show that these integrals 
can be expressed as quadratic combinations of the integrals H, A, and B. 
Therefore the following relations are valid: 

{A, C} p = aA 2 + f3B 2 + 2 7 AB + 8 A + eB + ( (23) 

and 

{B, C} P = aA 2 + bB 2 + IcAB + dA + eB + z (24) 

By taking an appropriate rotation of the integrals A and B we can always 
consider the case (3 = 0. 

The existence of the Poisson algebra (I2"3"l) - (j2"4"j) is not evident. The above 
form was assumed as obvious by several authors [I] -[7], [I7J I n Classical 
mechanics the Poisson algebra was not considered as an important point, but 
in Quantum Mechanics the existence of Poisson algebra permits the algebraic 
treatment of the superintegrable system [20]. In this paper we prove that the 
superintegrable systems can be classified using the properties of the Poisson 
algebra. The superintegrability is a global property of the system, and this 
fact is reflected in the Poisson algebra structure, which is indeed a global 
property. The mathematical proof of the existence of the algebra (|23p - (j2"4"j) for 
the two dimensional superintegrable systems can be found in the Appendix 

m 

The Jacobi equality for the Poisson brackets induces the relation 
{A {B, C} P } P = {B, {A, C} p } p 
The following relations 

b = —7, c = —a and e = —5 

must be satisfied. 

The integrals A, B and C satisfy the quadratic Poisson algebra: 

{A,B} P = C 

{A,C} P = aA 2 + 2 7 AB + 8 A + eB + ( (25) 
{B,C} P = aA 2 --fB 2 -2aAB + dA-5B + z 
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where a, 7, a are constants and 



z 



5 



c 

d 



e 



8(H) 

e{H) 

C{H) 

d(H) 

z(H) 



5 + 5^ 
eo + aH 



Co + Ci# + C2H 2 

do + d\H 

z + z 1 H + z 2 H 2 



where Si, e,, Q, di and Zi are constants. The associative algebra, whose the 
generators satisfy equations (12"5"I) . is the general form of the closed Poisson 
algebra of the integrals of superintegrable systems with integrals quadratic 
in momenta. 

The quadratic Poisson algebra (1251) possess a Casimir which is a function 
of momenta of degree 6 and it is given by: 



Therefore the integrals of motion of a superintegrable two dimensional system 
with quadratic integrals of motion close a constrained classical quadratic 
Poisson algebra (1251) . corresponding to a Casimir equal at most to a cubic 
function of the Hamiltonian (|26|) . 

In the general case of a superintegrable system the integrals are not neces- 
sarily quadratic functions of the momenta, but rather polynomial functions of 
the momenta. The case of the systems with a quadratic and a cubic integral 
of motion are studied by Tsiganov [31] [32] . 

IV Classes of superintegrable systems on a 
two dimensional manifold with quadratic 
integrals of motion 

The main result of the previous section is that the definition of the Casimir 
of the Poisson algebra, given by equation ff26l) determines uniquely the Pois- 
son algebra. This Poisson algebra is specific for each superintegrable system, 



K = C 2 - 2aA 2 B - 2 1 AB 2 - 26 AB- 
-eB 2 - 2CB + laA 3 + dA 2 + 2zA 
= k + k 1 H + k 2 H 2 + k 3 H 3 



(26) 



Obviously 



{K,A} p = {K,B} p = {K,C} P = 
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therefore it can be used for the classification of the possible superintegrable 
systems. Usually the proposed classifications of superintegrable systems as- 
sumed the definition of the manifold metric and the superintegrable systems 
were fixed for the given metric. In this paper we propose a classification 
which is based on the Poisson algebra. Let us consider a superintegrable 
system, which is described by a Hamiltonian H and two integrals of motion 
A and B. The integrability of the system imposes several choices which are 
determined two by Classes I and II of integrable systems, as it has been 
shown in section [Til These classes of super integrable systems are: 

Class I This class contains superintegrable systems, whose manifold metric 
and integrals of motion are written in a specific coordinate: 

ds 2 = g(£, rfidtdri, g{£, rf) = F(£ + rj) + G(£ - rf) (27) 
F(£ + r/) + G(£ - ??) + F(£ + ry) + G(£ - r?) 1 ] 



and 

A - t,2 i „2 _ O F(fj+ri)-G(Z-r,) , 

, a m+v) G(g-*>) - g(g-r?) 
^ + G«-i7) 

The second integral of motion has the general form: 



(29) 



b = a(OpI + B( v )pl - 2 p^ f ^ + ^_Ig^- v ) + g(e ' v) (30) 

where A(£) and B(rj) are not zero. 



Class II This class contains superintegrable systems, whose manifold metric 
and integrals of motion are written in a specific coordinate: 

ds 2 = <?(£, V )dtdri, <?(£, V ) = F( V )£ + G( V ) (31) 

ff = P&y f(v)Z + 9(v) ( o ?] 

F(77) C + G(77) + F( V )Z + G( V ) 1 ) 



and 



(33) 
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The second integral of motion has the general form: 




(34) 



where A(£) and B(rj) are not zero. 

The existence of a second integral of motion B implies that there is another 
Liouville coordinate system (X, Y) corresponding to the pair of integrals H 
and B. In this system there are analytic formulas for the second integral 
of motion by using the results of Proposition [TJ The analytic calculation of 
these formulas will be given in the following subsections. 

The Class I superintegral system is described by the integrals of motion 
which are given by equations (12TH3"Uj) . Correspondingly, the Class II super- 
integrable system is described by the integrals of motion which are given by 
equations fl3Tl - 134"|) . These integrals of motion should satisfy the relation 
(I26p . which is an identity relation between polynomials of sixth degree for 
the momenta. The coefficients of p| and p® in f|2oT) vanish and the following 
identities are true: 



by the solution of the above equations (1331) . 

The superintegrable systems on a manifold can be classified by the pos- 
sible solutions of equations (|35|) . The integral of motion A has a standard 
form, given by equation (|29|) . this form is the Liouville form and the coeffi- 
cients of pi and p 2 y are equal to 1. The second integral of motion B can be 
replaced by any combination of the form: 



where g, r, s are arbitrary constants. From this fact we can show that there 
are six subclasses of possible solutions of equation fl33|) . 




B 



qB + rH + sA 



• Subclass Ii This class corresponds to 



7 = 0, a = 0, a ^ 



if we choose a = — 6, then 
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Subclass I 2 This class corresponds to 

7 = 0, a^O, a = 
if we choose a = 8, then 

Subclass I3 This class corresponds to 

7 ^0, a = 0, a ^ 

if we choose 7 = 2, then we can show that all cases are equivalent to 
the choice 



A(0 = (e 5 + e-«) 2 , 5(7?) = (We"") 2 

and 

• Subclass IIx This class corresponds to 

7 = 0, a — 0, a = 

then 

A(0 = 1, S(t?) = 1 

• Subclass II2 This class corresponds to 

7 = 0, a = 0, a ^ 
if we choose a = — 6, then 

A(£)=£, 5(^=77 

• Subclass II 3 This class corresponds to 

7 = 0, a ^ 0, a = 
if we choose a = 8, then we can show that 

A(0=£ 2 , B( V )= V 2 

These classes of solutions will be studied in detail in Section [0 
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IV-a Class I superintegrable systems 

Starting from the definition of functions A(£) and B(rj) we can solve equations 
(IBT- fTTj) and the superintegrable system is fully determined. This procedure 
will be sketched in detail in the next paragraphs. 

Let us start by the known solutions B{rj) of equations (j551) . 

Equation (jSj) is written 

(A"(t)-B»( V ))(F(t + V )+G(t-ri)) + 

+3A'(fl (F'it + V )+ G\£ - v)) ~ 3fl'fa) (F'(£ + 17) - -»?))+ (36) 
+2 (A(0 - B( v )) (F"(Z + 77) + G"{£ -v)) = 

In the next paragraphs we will show that, the above equation can be sep- 
arated in two second order differential equations for the involved functions 
F(u) and G(v). The general solution of these equations are given by: 

F{u) = AiFi(u) + X 2 F 2 {u) + X 3 F 3 (u) + A 4 F 4 («) 
G(v)=£ 1 Gi{v) + £ 2 G 2 (v) + £ 3 G 3 (v) + UG^v) 1 ' 

where F^iu) and Gk{u) are functions which are not generally partial inde- 
pendent solutions of two second order differential equation (JHJ). Between the 
eight parameters A& and £k for k = 1, 2, 3, 4 only four among of them are 
independent. 

After the calculation of the functions F(u) and G(v) we can calculate the 
function (3(x, y) from equation The general form of the potential V(£, 77) 
is given by equation ffT4l) . 

After some elementary (but rather lengthy) algebraic calculation we can 
show that equation ffTTj) leads to a differential equation for the functions f(u) 
and g(y) which are involved in the definition (|l4p of the potential: 

fit + 77) {-3 B'( V ) (F'(£ + 77) - G'(Z - 77)) + 3 A'{$ (F'(£ + 77) + G'(£ - r,)) + 
+2 (A(£) - B( V )) (F"(£ + 77) + G»{£ - 

+9(Z - V) {-3 B'( V ) (F'(£ + T))- G>(£ - 77)) + 3 + + - 17)) + 

+2 (A(£) - 5(77)) (F"(£ + V) + G»(£ - 17))}- 

- (F(£ + 77) + G(£ - 17)) {-3 B'( v ) (/'(£ + 77) - - 77)) + 

+3 (/' (C + rj) + Sf{£ -r))) + 2 (A(0 - B(r])) (/"(£ + V ) + g»(t - 77))} = 

We can eliminate in the above equation the functions F(u) and G(v), which 
satisfy equation (l3"B"j) and finally the functions involved in the definition of 
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the potential satisfy the following equation: 

(A"(0-B"(v))(f(Z + v) + g(Z-v)) + 

+3A'(0 (g'(Z + V )+ g\{i - V )) - 3B'( V ) (<?'(£ + V ) - - r,)) + (38) 
+2 - B(r])) (/"(£ + »7) + <f - i?)) = 

This equation is indeed the same as (1361) and the general solution has been 
given by equation (|37|) . Then we have proved the following general proposi- 
tion: 

Proposition 2 The general form of the potential V(£, g) in Liouville coordi- 
nates for a superintegrable system of Class I is given by the general formula: 



F(i + ,) + G({ - 7,) 



where both the pairs of functions f(u), g{v) and F(u), G{v) satisfy the same 
differential equation /[3b]) and §E3j) - 

Therefore the solutions of equation (l38l are given by: 

f(u) = piFi(u) + p2F 2 (u) + p 3 F 3 (u) + piF 4 (u) , , 

g(v)=r l G 1 (v) + r 2 G 2 (v) + r 3 G 3 (v) + r 4 G A (v) ^ 

The equations of motion don't depend on the shifts of the potential by one 
constant, this fact explains physically the identity of differential equations 
(1361) and (EE}. 

The calculation of the second integral of motion B is now straightforward. 
If the functions A(x) and B(y) are fixed as solutions of the characteristic 
equations (1351) then we can introduce the new coordinates: 

X= [*L= and Y= (40) 



\[Mx) ' J yfm) 

In these coordinates the metric of the manifold can be written 

ds 2 = g(x,y)dxdy = g(X,Y)dXdY (41) 

where 



g(X,Y)=g(x,y) yj 1 A{x) B{y) 
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and 



H = —, 7 + V{x, y) = + V(X, Y) 



The integrability of the above Hamiltonian implies: 

g(X, Y) = F(X + Y) + G(X - Y) (42) 

The above relation implies: 

g(^f,^) = F(U)+G(V) 
g(u±d^u^ = F (u)+G(d) 

sfe^) = F(c) + G(V) 



9(f,^)' = F(c) + G(d) 

where c and d are two arbitrary constants. Therefore the functions F(U) 
and G(V) are calculated up to one constant by: 

m=9(^^)-\g(^ c -^) + » ,,,, 

where \i can be an arbitrary chosen constant. Starting from the well known 
form of the potential: 

V(x,y) = W X,V , w(x,y) = f(x + y) + g(x - y) 

9{x,y) 

we can show that 

V(x,y) = V(X,Y) = ^ll 

_ 9{X,YY (44) 

^(X, Y) = f(X + Y)+ g(X — Y)= w(x, y) J A(x) B(y) 



The functions f(U) and g(V) are calculated by: 
Then the second integral of motion in the coordinates X, Y is: 

^ + 2PXPy F(X + Y)+G(X-Y) + 

1 /(x+y)G(x-y)-g(x-y)F(x+y) 
F(x+y)+G(x-y) 

After calculating the above integral in the coordinates X, Y, we can compute 
analytically the second integral B in the original coordinates x, y. 



(45) 
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IV-b Class II superintegrable systems 

Equation (jSJ) is written 

(A"(x)-B"(y))(F(y)x+G(y)) + 

+3A'(x)F(y) - 3B'(y) (F'(y)x + G'(y)) + (46) 
+2 (A(x) - B{y)) (F"(y)x + G"(y)) = 



(47) 



The general solutions of equation fl46|) are given by: 

F(y) = \ 1 F 1 (y) + \ 2 F 2 {y) + A 3 F 3 (y) + X 4 F 4 {y) 
G(y)=£ 1 G 1 (y) + £ 2 G 2 (y) + £ 3 G 3 (y) + £ 4 G 4 (y) 



where F k (y) and G^{y) are partial independent solutions of two second order 
differential equations with several constant parameters. Between the eight 
parameters and £k for k = 1, 2, 3, 4 only four among them are linearly 
independent. 

After the calculation of the functions F(y) and G(y) we can calculate the 
function fi(x, y) from equation (TT5|) . The general form of the potential V(x, y) 
is given by equation 0191) . After some elementary (but rather complicated) 
algebraic calculation we can show that, equation (120]) leads to an differential 
equation for the functions f(y) and g(y), which are involved in the definition 
( THE]) of the potential: 

{A"{x)-B"{y)){f{y)x + g{y)) + 

+3A'(x)f(y) - 3B'(y) (f'(y)x + g'(y)) + (48) 
+2 (A(x)-B(y))(f"(y)x + g"(y))=0 

Then we have proved the following general proposition: 

Proposition 3 The general form of the potential V(£, rj) in Liouville coordi- 
nates for a superintegrable system of Class II is given by the general formula: 



F(t,)i + G(ri) 



where both the pairs of functions f(u), g(v) and F(u), G(v) satisfy the same 
differential equation fijh] ) and 



Therefore the solutions of equation (1481) are given by: 

f(y) = p 1 F 1 (y) + p 2 F 2 (y) + p 3 F 3 (y) + p 4 F 4 (y) 
g(y) = nG x {y) + r 2 G 2 (y) + r 3 G 3 (y) + r 4 G 4 (y) 
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The equations of motion don't depend on the shifts of the potential by one 
constant, this fact explains physically the identity of differential equations 
([H (or flUD ) and (J3§D (or (PJ). 

The calculation of the second integral of motion B is now straightfor- 
ward, we can use the same procedure of solution as it has been described by 
equations fjlOl - HoD 

V Classification of two dimensional superin- 
tegrable systems with two quadratic inte- 
grals of motion 

In this section we give the analytical solutions for the different classes of 
superintegrable systems. As we have shown there are two general classes of 
superintegrable systems, each class has 3 subclasses. 

V-a Class I superintegrable systems 
V-a.l Subclass Ii of superintegrable systems 



A(0 = f , B{rj) = V 



F(u) = AXu 2 + ku + i//2, G(v) 
f(u) = A£u 2 + ku + n/2, g{v) 



-Xv 2 + fi/v 2 + u/2 
-£v 2 + m/v 2 + n/2 



(50) 



ds 2 = g(£, n) dr), g(£, rj) = F(£ + V ) + G(£ - V ) 

H=-%pT + V{£,T,) V(^ V ) = ^l, w(£,T,) = f{Z + T,)+g(t-r,) 
The other integral of motion is: 

F(Z + v)-G(£-ti) 



A=p^+p r/ - 2ptf v 



J(Z + v)G(t-y) - 9(Z-V)F(Z + V ) 
F(C + V) + G^- V ) 
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We introduce the functions: 



~, . Am 6 k,u 4 vu 2 a 

G(t,) = -^--^--^ + 4 
v ; 256 128 16 v 2 

~, . lu & ku 4 nu 2 m 

i[u) = ^56" + ~m + ~lo~ ~ 

. . Zt> 6 /cv 4 no 2 m 

^ " _ ~256" ~ ~1~28 16~ + 

The second integral of motion is: 

yx _ ^ y _ F(X + y) + G(X - Y) 
n f(X + Y)G(X - Y) - g{X - Y)F(X + Y) 
F(X + y) + G(X — Y) 

where 

The constants of the Poisson algebra are: 

a = = 7 = 5 — 16(kH - k) e = 256(A# - t) 
C = -32(kH - k)(uH - n) a = -6 d = S(uH - n) 
z = 8(uH - n) 2 - 128(A# - £)(jiH - m) 

K = 32(uH - n) 3 + 512(XH - f)(jiH - m)(vH - n) - Q4(kH - k) 2 (fiH 
V-a.2 Subclass I 2 of superintegrable systems 



, ■ ) k n , o / / ' 7i 
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ds 2 = gfo V ) di d V , gfo rj) = F(£ + V ) + G(£ - rj) 

H = ji~n) +V ^ 7]) m " ) = W$' v&i) = m + n) + 9(t-n) 

The other integral of motion is: 

A 2 M 2 9 F(Z + v)-G(Z-ri) 

J{Z + v)G(Z-y) - g(Z-v)F(£ + ri) 

F(Z + V ) + G(Z- V ) 

We introduce the functions: 



F(u) = 4 A e 2u + v e M , G(v) = - 6 2 + ' 



V 



f(u) = 4£e 2u + ne u , g(v) = — ^ -t 



(l + e v j z ' (-1 + e^) 2 



The second integral of motion is: 



(53) 



x _ y F(x + y) + G(x - y) 

| 4 /(X + y)G(X - Y) - g(X - Y)F(X + Y) 
F(X + Y) + G(X - y) 

where 

X = ln£, Px = tPt, Y = \nr], p Y =VP v 
The constants of the Poisson algebra are: 

a = 8 /3 = 7 = 5 = e = 256(A# - t) 

( = -32(uH - n) 2 + 256(A# - i) ((// - re) if - (m - fc)) 

a = d = z = 32 ((k + -(fc + m)) (v H - n) 

K = 256 (A #-£) ((k + //)// - (k + m)) 2 +128 ((« - /x)# - (A; - m)) (z/# - n) ; 
V-a.3 Subclass I 3 of superintegrable systems 

A(0 = (e« + e"«) 2 , 5(^) = (e" + e^) 2 
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„, x ne 2u \e u (l + e 2u ) . . ue 2v v e v (1 + e 2v ) 
F(u) = ~ H G(v) = , H ^ 

v ' / 1 , „9->,\2 / 1 i _9«\'i ' v > I 1 , _9-,i\2 / 1 , _9i>\2 



(-l + e 2 «)" (-l + e 2 ")' v/ {-l + e 2 -y (-1 + e 2v Y 

ke 2u £e u (l + e 2u ) me 2 ^ ne^fl + e 2 ^) 
/(u) = * H i o(?j) = ^ H i ^ 

(54) 

cfe 2 = 77) di d V , g(£, 77) = F(£ + 77) + G(£ - 77) 
The other integral of motion is: 



1 /'i • 'V Jyr- — ; - : 

• 4 



+ 77) + G(f - 77) 
/(£ + 77)^-77) - ff(e-r7)F(e + ?7) 



(55) 



F{i + V) + G(e - 77) 

We introduce the functions: 

~. . (k + 2 A) 9 2u — a 9 A + z/ 

F(u) = tan 2 u + — - cot 2 u + 

4 4 2 

G(?j) = - - tan 2 v + — cot 2 v H 

4 4 2 

~, , (k + 2£) o 2n -m 9 £ + n 

f(u) = - — - — - tan 2 m + cot 2 u + 

4 4 2 

N (2£-k) 2 m + 2n 2 ^ + ™ 

q[v) = tan v H cot v H 

1 ' 4 4 2 

The second integral of motion is: 

2 2 f(x + r) - G(x - r) 

B= p x +p Y -2p x p Y ^ '- J '-+ 

X _ y F(X + y) + G(X z Y) 

| 4 /(X + Y]G{X - Y) - g(X - Y)F{X + Y) 
F(X + Y) + G(X - y) 

where 

X = arctan(e^), = (e^ + e~^) p^, Y = arctan(e ,? ), Py = + e~ v ) p v 
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The constants of the Poisson algebra are: 

a = -32 = 7 = 8 5 = e = ( = -32 (XH - £) (vH - n) 

a = d = 64 (k - m) - 64 (k - n) H 

z = 32 ((A - u)H -{£- n) f - 32 {{kH - k) (/iH - m) 

K = Q4 (kH - k) (uH - nf - QA(XH - tf (fiH - m) 

V-b Class II superintegrable systems 
V-b.l Subclass IIx of superintegrable systems 

A(0 = i, b( v ) = i 

F(?7) = kt] + X, G(ri)=nri + u ^ 
f(r)) = kr] + £, g(r]) = mr] + n 

ds 2 = g(£, 77) di d V , g(£, r,) = £ Ffa) + 
The other integral of motion is: 

A = _ 2n P j f M in 2 g /(„) + g („)) j £W *, + 2 / /( ^ 

We introduce the functions: 



g ( „ ) = _^ + (A^), + , 

-/ N /cm 2 + m) -u n 
/(M) = ^ + ^^ + 2 

. . fcw 2 (£ — m) v n 



(57) 
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The second integral of motion is: 



B= pl + p 2 v -2pzp v 



F(Z + V)-G(S-V) 



+ 



+ 4 



f(Z + v)G(Z-v)-gjZ-y)F(Z + Ti) 



The constants of the Poisson algebra are: 

a = p = 7 = 6 = 8(k-KH) e = ( = 8 (\H 
a = d=16(k-KH) z = 8 (\H - if - {fiH - mf 

K —16 (uH - nf {kH - k) - 32 (\H - i) (fiH - m) (uH 
V-b.2 Subclass II 2 of superintegrable systems 



A(0 = t, B(n) = v 

F(r]) = — + X, G{r]) = 3k^+ Xr] + ^- + v 
ds 2 = gfor,) did V) g&n) =ZF(n) + G(rj) 



if 
n) 

(58) 




2PtPr,f F( V ) dr) 2 (g f(rj) + g(rj)) J F(rj) d V 

g(t,,v) g(t,,v) 



+ 2 
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We introduce the functions: 



~, . Am 4 ku 3 vu 2 liu 
F ( u ) = T^T + ~T7T + —7T + 



G(v) 



128 16 16 4 

Af 4 KV 3 fiV vv 2 

"128 ~16~ ~4 16~ 



~. . £m 4 /cm 3 nu 2 mu 



128 16 16 4 

f 4 kv 3 mv nu 2 



= + -T7T- + 



128 16 4 16 
The second integral of motion is: 

^ X _ ^ y _ F(X + Y) + G(X - Y) 
n f(X + Y)G(X -Y)- gjX - Y)F{X + Y) 
F(X + Y) + G(X - Y) 

where 



x = 2 y£, Px = y/ZPt, Y = 2^/rj, p Y = VvPv 
The constants of the Poisson algebra are: 

a = (3 = 7 = 5 = 4(£-XH) e = ( = 8(kH - k) 2 
a = -6 d=8(uH-n) z = —8 (kH — k) (liH — m) — 2 (vH - 

K = 8 (A H - £) (fiH-m) 2 -16 (k H - k)(n H - m) (u H - n) 
V-b.3 Subclass II 3 of superintegrable systems 



m = e, b( V ) = v 2 

F( V ) = \ V + - 3 , G( V ) = u+-^ 
f(rj) = £r] + g(r])=n + ^ 



29 



ds 2 = gfor,) d£d V , g(^rj) =£F( V ) + G(rj) 

The other integral of motion is: 

A _ p 2 2p i p,JF(r ] )dr ] 2 (U( V ) + g( V )) j F( V ) d V | g r 



We introduce the functions: 



F(u) = Xe 2u + ue u , G(y) = re e 2v + fi e v 
f(u) = £e 2u + ne u , g(v) = ke 2v + me v 



(61) 



The second integral of motion is: 



, n F(X + Y)- G(X - Y) 

B= p\ + p\ - 2p x p Y '- J 

X _ ¥ F(X + Y) + G{X - Y) 

| ^{X + Y)G{X -Y)- g{X - Y)F{X + Y) 
F{X + Y) + G{X - Y) 

where 

X = ln£, Px = iPi, Y = \nr], p Y =VP v 
The constants of the Poisson algebra are: 

a = 8, = 0, 7 = 0, 5 = 0, e = 0, ( = 32 (kH — k) (XH — £) , 
a = 0, d = 0, z = 32 {\iH -m){yH - n) 

K = 64 (A H - t) (fiH - mf - 64 (k if - Jfe) (uH - nf 

All the above superintegrable systems generally are defined on manifolds 
which have neither constant curvature nor are they surfaces of revolution. 
All the known superintegrable systems are defined on manifolds of constant 
curvature or on surfaces of revolution. Therefore we have proved that there 
are new superintegrable systems, which have not yet been studied. 
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VI Superintegrable systems corresponding to 
Koenigs essential forms 



Class Ii Using the coordinate transformation 

1 



X 2 , Pt: 



Pa 



X 



V 



y 2 , p v 



Py 

y 



the metric of the Class I 2 superintegrable systems is reduced to the 
metric of the essential form VII. 4 [IJ vol IV, p. 385], if 



K 



16A, 



A = 16A, 



-An 



V 



AA 1 



The corresponding superintegrable system (using the coordinates of 
Ref. [1]) is given by the Hamiltonian: 



H 



p x p y + w(x,y) 

g(x,y) 



g{x,y) = A Q 
+ A 2 

w(x,y)= a 
+ a 2 



+ A X [(x + yf - (x - yf] + 
(x + yf - (x - yf] + A 3 [(x + yf - {x - yf] 



i 

[x+y) 2 



(x-y) 2 



{x+y? j^yYj + ai [(x + yf-(x-yf] + 
x + yf -{x- yf] + a 3 [(x + yf - (x - yf] 



where only three of the constants ao, aj, a 2 , a 3 are independent, i.e. we 
can put one among them equal to zero. Using relations (51) we have 
that in x, y coordinates the other integral of motion is: 



A(x,y) 



2 Px 2 Py 



+ 



„ „ ®(x,y)-$(x,y) 
PxP y <S>(x,y)+^(x,y) 
2 Mx, y)^(x, y)-tp(x, y)$(x, y) 
^(x,y)+^(x,y) 



where 



$(x,y) 
*(x,y) 
(f)(x,y) = 
tp(x,y) 



A 



i 



(x+y) 2 

A 1 



a 



(x-y) 2 
1 



+ Ai(x + yf + A 2 (x + yf + A 3 (x + yf 
- yf - A 2 (x - yf - A 3 (x - yf 



{x+y ) 2 + a ^ x + y) + a ^ x + v) + a s( x + y) 

1 a x {x - yf - a 2 (x - yf - a 3 (x - yf 



-a 



(x-y) 2 
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while using relations (50) we have that in x, y coordinates the second 
integral of motion is: 



B(x,y) 



— 1 „2 I 12 2 w " $ ^.!<)-*(^i/) 
X 2F X ~ y 2 ry xy ^(x,y) + ^(x,y) 

, ^ </>(a:, y)^(x, y)~'tp(x, y)$(x, y) 
x y Q(x,y)+$!(x,y) 



where 



+ A 3 ( 



(x + yf - {x-y)' 
(x + yf -{x- yf 



+ A 2 (x + yf - (x - yf 



+ 



+ (x + yf(x-yf- (x + yf(x-yf) 



V{x, y) = A 
<j>(x, y) 



+ ~A X (x + y) 2 -(x-yY 



(x+y) 2 (x-y) 2 _ 

- A 3 [(x + yf(x - yf - (x + y) 2 (x - yf) 



{x + yf - (x-y)' 
(x + yf -(x- yf 



+ a 3 ( 

+ (x + yf(x-yf- (x + yf(x-yf) 



+ a 2 



(x + yf - (x - yf 



+ 



if>(x, y) = a 



+ \ai (x + yf-(x-yf 



(x+y) 2 (x-y) 2 _ 

«3 {(x + yf(x - yf -(x + yf(x - yf) 



Class I2 Using the coordinate transformation 
1 



f = __ CO s(2a;), p^ = . Px , 
2 sm(2x) 



1 p v 
r} = --coa{2y), p v = — — , 

2 sm(2y) 



the metric of the Class I2 superintegrable systems is reduced to the 
metric of the essential form VII. 4 P, vol IV, p. 385], if 



K = Ai, X = — 8A 3 , fi = —A , v 



-2Ai 



The corresponding superintegrable system (using the coordinates of ref 
p]|) is given by the Hamiltonian: 



H 



p x p y + w(x,y) 
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where 
9(x,y) 



= A 
+ A 2 



sin 2 (x + y) sin 2 (a; — y) 
cos 2(x + y) — cos 2{x — y) 



w(x,y) = a 
+ a 2 



sm 2 (x + y) sin 2 (x — y) 
cos 2{x + y) — cos 2(x — y) 



' 1 cos 2 (x + y) cos 2 (x — y) 
+ A 3 [cos4(x + y) — cos4(x — y)} 



cos 2 (x + y) cos 2 (x — y) 
+ a 3 [cos4(x + y) — cos4(x — y)] 



+ 



+ 



where only three of the constants ao, a%, a 2l a 3 are independent, i.e. we 
can put one among them equal to zero. Using relations (1531 we have 
that in x, y coordinates the other integral of motion is: 



A(x, y) ■- 4 
+ 



\cot 2 (2x)p 2 x + \cot 2 (2y)p 2 y - 
i 



PxPy 



$(x,y)-^(x,y) 



4"" \"i>lfy 2tan(2x)tan(2y) <&(x,y)+9(x,y) 
x, y)^(x, y)-ip(x, y)<b(x, y) 



tan(2x) tan(2y) 



<S>(x,y)+^(x,y) 



where 



y) = A 2 [cos2(x + y) — cos2(x — y)\ + A 3 [cosA(x + y) — cos4(x — y)} 



V(x,y) = A 



i 



i 



sin 2 (x+y) sin 2 (x—y) 



A, 



i 



i 



cos 2 (x+y) cos 2 (x—y) 



(x,y) = a 2 [cos2(x + y) — cos2{x — y)] + a 3 [cosA(x + y) — cos4(x — y)] 



$(x,y) = a 



+ ai 



i 



i 



cos 2 (x+y) cos 2 (x—y) 



sin 2 (x+y) sin 2 (x—y) _ 

while using relations fl52|) we have that in x, y coordinates the second 
integral of motion is: 



B(x, y) 



1 —v 2 + 

2x)t ) x T 



where 
V(x,y) 



+ 4 



At 



<5>(x,y)-f(x,y) 



sin 2 (2x)fx ' sin 2 (2y)fy 2 sin(2x) sin(2y) &(x,y)+ty(x,y) 

1 <t>(x, V) ^(x, y)-ip(x, y) y) 



sin(2x) sin(2y) 



1 



<3>(x,y)+V(x,y) 



1 



cos 2 (x+y) cos 2 (x—y) 



+ \A 2 [cos2{x + y) — cos2{x — y)} + 



+ 4A 3 cos 2 (x + y) cos 2 (x — y) [cos2(x + y) — cos2{x — y)} 



A, 



i 



i 



sin 2 (x+y) sin 2 (x—y) 



+ \A 2 [cos2(x + y) — cos2(x — y)] — 



4 A 3 sin 2 (x + y) sin 2 (x — y) [cos2(x + y) — cos2(x — y)] 
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<f>(x,y) 



a i 



cos 2 (x+y) cos 2 (x—y) 



+ \a 2 [cos2(x + y) — cos2(x — y)] + 



+ 4a 3 cos 2 (x + y) cos 2 (x — y) [cos2{x + y) — cos2{x — y)} 



ip{x,y) 



a 



i 



sin 2 (x+y) sin 2 (x—y) 



+ \a 2 [cos2(x + y) — cos2(x — y)] — 



Class Is 



— 4 a 3 sin 2 (x + y) sin 2 (x — y) [cos2(x + y) — cos2(x — y)] 
Using the coordinate transformation 



e = in 

7] = In 
where 



A 

{x)~p{u) X ) 



p{y)-p{u\) 
A 

p{xy)-p(ui) 



+ -ln 

2 



P{u\)-p{u2) 

A 

p(^i)-p(^ 2 ) 



+ 1 



+ 1 



In 



p{ui)-p(u2) 

A 

P{ui)-P{u2) 



+ 1 



Pi, 



p{2x) - pjujt) 
2A 

2A 



Pi 



A^ 



- K w 2))(p(wi) - p(o; 3 )) 



the metric of the Class I3 sup erintegr able systems is reduced to the 
metric of the essential form VII. 1 P, vol IV, p. 385], if 



k = 2(A 2 + A 3 ), \ = A 2 -A 3 , fi = -2(A + A 1 



-Aq + A x 



The corresponding superintegrable system (using the coordinates of ref 
[I]) is given by the Hamiltonian: 



H 



PxP y + w(x,y) 

g(x,y) 



where 



g(x, y) = A (p(x + y) - p{x - y)) + A x (p(x + y + to x ) - p(x - y + 
+ A 2 (p(x + y + u 2 ) -p(x-y + u 2 )) 
+ A 3 (p(x + y + u 3 ) -p(x-y + u 3 )) 

w(x, y) = a (p(x +y)- p(x - y)) + a x (p(x + y + U) X ) - p(x - y + wi)) 
+ a 2 (p(x + y + lu 2 ) - p(x - y + uj 2 )) 
+ a 3 (p(x + y + uj 3 ) - p(x-y + u 3 )) 

where only three of the constants a , a±, a 2 , a 3 are independent, i.e. we 
can put one among them equal to zero. Using relations fl54l) we have 
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that in x, y coordinates the other integral of motion is: 



A(x, y) = £s(p(2x) - pM)pj + 4sr(p(2j/) - pM)p$- 



y) = A 2 0(x + y + ^ 2 ) - p(x-y + u 2 )) + 
+ A 3 (p(x + y + u 3 ) - p(x - y + u 3 )) 

V(x,y)= A (p(x + y) - p(x - y)) + 

+ Ax (p(x + y + u x ) - p{x - y + u x )) 

(f)(x,y)= a 2 (p(x + y + uj 2 ) - p(x - y + u 2 )) + 
+ a 3 {p{x + y + u 3 ) - p(x-y + u 3 )) 

ip(x,y)= a (p(x + y) - p(x - y)) + 

+ ax (p(x + y + ui) - p(x - y + ui)) 

while using relations (|55l) we have that in x, y coordinates the second 
integral of motion is: 



+ 




where 



B(x, y) 



+ 



1 (,„(n„.\ , 




+ 




(p(w8)-p(wi)) (p(w2)-p(o) 3 )) *(a;, y)+9(x, y) 



where 



y) 



A a + y + wi) - p(x - y + wi)) + 
^2(^(2 + 2/ + w 2 ) - p(x - y + w 2 ))- 



+ 




*(ar, y) 



A (p(x + y)-p(a;-y)) + 

^(pfa + y + w 3 ) - p(ac - y + w 3 ))+ 



+ 



+ 



i (A -4- A -1-4 -i- A \ V P(2a:)-p(^3) V p(2j/)-p(^ 3 ) 
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$(x,y)= ai (p(x + y + - p(x - y + ui)) + 
+ a 2 (p(x + y + u 2 ) - p{x - y + u 2 ))- 



1 / L „ j i \ V p(2x)-p(^ 3 ) v/ p(2j/)-p(g;3) 

ip(x,y)= a (p(x + y) - p(x - y)) + 

+ a 3 (p(x + y + lu 3 ) - p(x - y + uj 3 )) + 



If i i _|_ n ^ V p(2a:)-p(^3) V p(2j/)-p(^ 3 ) 

2 + Blt02 + (p(^ 3 )-p(^l))(p(a;2)-p(a; 3 )) 



Class Hi This case is not covered by Table VII of Koenigs pQ, Vol. IV, p. 
385J. This class corresponds to the Kress [T6] equivalence class [0, 11] of 
the nondegenerate superintegrable systems En, E 20 of ref [6]. 

The Hamiltonian is 

_ PtPri + k^T] + + mrj + n 
n£r] + A£ + fir] + v 

the integrals of motion are 



(f rf + A77) (^77 + ^ + 

mrj + 



+ A£ + fir] + f 
k^t) + A£ + fir] + f 



2 



f (£ 2 + + A77 + /i£J (fcf 77 + ^ + mrj + n 
k£?7 + A£ + pir] + f 



The case VI 6 in Table VI[H Vol. IV, p. 384] Koenigs studied separately 
the cases where k = and k^O, A = /x = 0. 



Class II2 Using the coordinate transformation 



> 1 2 1 2 Pj/ 

the metric of the Class II 2 superintegrable systems is reduced to the 
metric of the essential form VII. 5 [U vol IV, p. 385], if 

k = 2y/2A u A = I6A0, At = V2A 3 , v = 4A 2 
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The corresponding superintegrable system (using the coordinates of ref 
[Tj) is given by the Hamiltonian: 

= p x p y + w(x,y) 

where 

g(x, y) = A [(x + yf - (x - yf] + A x [{x + yf - {x - yf] + 
+A 2 [(x + yf -(x- yf] + A 3 [(x + y) - (x - y)\ 



w(x, y) = a [(x + yf - (x - yf] + a x [(x + y f - (x - y) 3 ] + 
+a 2 [(x + yf - (x - yf] + a 3 [(x +y)-(x- y)) 

where only three of the constants a , a x , a 2 , a 3 are independent, i.e. we 
can put one among them equal to zero. Using relations ( l59l) we have 
that in x, y coordinates the other integral of motion is: 

Afrr nl \ — lr> 2 -I- In 2 — r> r> ®( x ,y)-y( x >v) _r_ o 4>(x,y)y(x,y)-4>$(x,y) 
Sly^, y) ^ 2 Py V*VV §(x,y)+9(x,y) ^ $(x,y)+$(x,y) 

where 

<£(x, y) = A (x + yf + A x (x + yf + A 2 (x + yf + A 3 (x + y) 

f (x, y) = -A (x - yf - A x (x - yf - A 2 (x - yf - A 3 (x - y) 

4>(x, y) = a (x + yf + a x (x + yf + a 2 (x + yf + a 3 (x + y) 
ip(x, y) = -a (x - yf - a x (x - yf - a 2 (x - yf - a 3 (x - y) 

while using relations (IB^|) we have that in x, y coordinates the second 
integral of motion is: 

B(x, y)= ®{x)dx -2$$ J <$>(x)dx + J <f>(x)dx 

where 

$(x) = lQA Q x + AA X 
(f)(x) = 16 clqx + 4 ai 
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Class II 3 Using the coordinate transformation 



t, = -\cos{2y), n = -^Pv, V = e**, P v = --^Px 

the metric of the Class II 3 superintegrable systems is reduced to the 
metric of the essential form VII. 3 P, vol IV, p. 385], if 

k = 2(A 1 +iA Q ), X = 2(iA -A 1 ), n = - - (A 3 +i A 2 ) , i/=-(A 3 -iA 2 ) 
where 

g(x, y) = A [sin4(x + y) — sin4(x — y)] + A\ [cos4(x + y) — cos4(x — y)] + 
+A 2 [sin 2{x + y) — sin 2{x — y)} + A 3 [cos 2{x + y) — sin 2(x — y)] 

w(x, y) = a [sin 4(x + y) — sin 4(x — y)\ + a x [cos A(x + y) — cos 4(x — y)} + 
+a 2 [sin 2(x + y) — sin 2(x — y)} + a 3 [cos 2(x + y) — sin 2(x — y)] 

where only three of the constants a , ai, a 2 , a 3 are independent, i.e. we 
can put one among them equal to zero. Using relations flBTl) we have 
that in x, y coordinates the other integral of motion is: 

A(x, V) = ~\Px + \ C0t2 {^y)p 2 y + 2 i tan(2y) ${ x [ ~~ 
_ 1 4>{x,y)y{x,y)-$(x,y)®(x,y) 

2 i tan(2y) i>(x, y)+^{x, y) 

where 

f>(x, y) = (A + i A) e 4i x sm(4y) + (A 2 + i A 3 ) e 2i x sin(2y) 
V(x, j/) = (A - i Ai) e~ 4i x sm(4y) + (A 2 - i A3) e" 2i x sin{2y) 

4>(x, y) = (a + i ai) e 41 x sin(4:y) + (a 2 + i a 3 ) e 21 x sin{2y) 
ip(x, y) = (a — i ax) e~ 41 x sin(Ay) + (a 2 — i a 3 ) e~ 21 x sin{2y) 

while using relations fl60l) we have that in x, y coordinates the second 
integral of motion is: 

where 

y) = -8A cos(Ax) + 8A x sm(4x) 
y) = — 8a cos(4x) + 8aism(4x) 
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VII Superintegrable potentials on a surface 
of revolution with two quadratic inte- 
grals of motion 

A manifold which is described by a metric of the form 

ds 2 = g{x + y)dxdy or ds 2 = g(x — y)dxdy 

is called a surface of revolution. 

The above condition is possible only for a specific choice of the parameters 
K, A, /i and v. In many cases the superintegrable systems can be calculated 
by using the general forms which are studied in Section |Vj The general forms 
of these systems by revolution in many instances are given by the formulas: 

P i p ri + /(£ + v) + g(Ji - v) PzP v + f(£ + v) + 9(£-y) 

F(C + rj) ° r G(£ - rj) 

and 

„ P xPY + f(X + Y)+g(X-Y) p xPY + f(X + Y)+g(X-Y) 

tl = — or n = — 

F(X + Y) G(X - Y) 

But we must notice that the Liouville or the Lie coordinates are not always 
the appropriate ones for concluding whether a surface is a surface by revo- 
lution. Among the parameters k, A, // and u the surfaces of revolution are 
determined by two independent parameters. In Table [2] these special values 
of the parameters k, A, fi and v are shown. We notice the corresponding 
potentials given in references [13] and [12] . This classification scheme shows 
that there is the case Ru which is not given in the above references and it 
is a new not known superintegrable system. A detailed description of the 
superintegrable systems on surfaces of revolution are given further in the 
present section. 



Rx: Class I± « = 0, A = 



ordinates is given by 
H 



The form of the Hamiltonian in Liouville co- 



(£-v) (£-»?) 
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Class 


K 


A 






Essential 
from [1] 


Classes 
from [16] 


h 


16A 2 


I6A3 


-A 




VILA 


[3,2] 


h 


Ax 


-8A 3 


-A 


-2A 2 


VII. 2 


[21,2] 


h 


2{A 2 + A3) 


A 2 -A 3 


-2(A + Ai) 


-A + A t 


VIII 


[111,1] 


Ih 












[0,11] 


Ih 


2^2A X 


16A 


V2A 3 


4A 2 


VII. 5 


[3,11] 


Ih 


2{A 1 + iA ) 


2(iA -A 1 ) 


-±(A 3 + iA 2 ) 


L 2 (A 3 -iA 2 ) 


VII.3 


[21,0] 



Table 1: Essential forms of Table VII in Ref. [T] and equivalence classes of ref. 















Potentials 


Potentials 




Class 


K 


X 




V 


by revolution 
from ref [13] 


by revolution 
from ref [12] 


Ri 


h 












2[A] 




R 2 
















(1) 


R 3 


h 












2[B] 




R4 














m 




R 5 














4[A] 




Re 




-fi 









2[C] 




R 7 


h 












4[B) 




R 8 




-// 


V 






4[C] 




Rg 


Ih 













(2) 


Rio 














m 




R11 


Ih 














new 


R12 


Ih 




— K 


—V 




3[D] 




Rl3 




-n 










m 





Table 2: Potentials by revolution with two quadratic integrals of motion 
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By the coordinate transformation 

v + iu v — iu 
f = — 2 — ' 11 = — 2 — >PZ = P V ~ l P u >Pv=P v + l Pu 

and putting fi = — 1, v — 1, the Hamiltonian 2 [A] of reference [13] is 
obtained 

The form of the Hamiltonian in Liouville co- 



R 2 : Class Ii A = 0, y, = 

ordinates is given by 

K (v + + ^ « + + f 

By the coordinate transformation 

£ = u + iv, r] = u-iv,p^ = ^{p u - ip v ) ,p v = -(pu + ip v ) 

and putting k = 1/2, v — the Hamiltonian (1) of reference [12] is 
obtained 

E = vl±vj , 16£(4 M 2 + ^) | n f | 

4ii Au u uv 2 



R 3 : Class I 2 k = 0, A = 0, 



dinates is given by 



The form of the Hamiltonian in Liouville coor- 



By the coordinate transformation 

v + iu v — iu 
£ = — ^ — ' V = — 2 — ' Pe = Pv ~ iPu 'Pv = P v + l P u 

and putting fj, = —1, u = 1, the Hamiltonian 2[B] of reference p3] is 
obtained 
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i?4 : Class I 2 k = 0, /i = The form of the Hamiltonian in coordinates 



(X, Y) is given by 

Px p Y + Ue 2 ^ + ne^ + k fJ Y + m ' " v ; 1 



4Ae 2(x+y) + ^(x+y) 
By the coordinate transformation 

1 / 4 \ 1/ 4\ 
X = In n \,Y= In 



2 \(iv + u) 2 J , 2 y (v + iu) )' 

Px = 2 + u ^ Pu + ( v ~ iu )Pv) > = ~~ ~~ ( u + 

and putting A = 1, z/ = 4 the Hamiltonian 3[S] of reference [13] is 
obtained 

A + u 2 + v 2 



R 5 : Class I 2 A = 0, v = The form of the Hamiltonian in coordinates (X, Y) 



is given by 

p x p Y + Ue 2 ^ + ne^ + k + m ' v ! 5 



,_, (i +e x+y)2 -r ( _ 1+e x+y)2 

K (i +e x+Yy + V (_ 1+e x+y)2 

By the coordinate transformation 

X = In Q(x - iy)^ , y = In Q(x + iy) 

Px = ^((x-i y)p x + (i x + y)p y ) ,py = ~((x + i y)p x + {y-i x)p y ) 

Putting k = ^j^, \i = -I 2 , the Hamiltonian 4 [A] of reference [13] is 
obtained 

4a: 2 !/ 2 (gj + pi + n + \ ((2 + q)fc + (a - 2)m) + 4,) + I (x 2 + y 2 )) 
H ~ y 2 (a-2) + x 2 (2 + a) + + ™ 
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ordinates is given by 



Rq\ Class I2 re = —fi, v = The form of the Hamiltonian in Liouville co- 



PvP( 



-(^ + A(r/ + 2 -A(e-r/) 2 + ^ 



J ^ + A(r7 + £) 2 -A(£-7 7 ) 2 



By the coordinate transformation 

1 1 
£ = u + iv,ri = u-iv,p ( : = -{p u ~ ip v ), Pr, = ^(Pu + ip v ) 

and putting A = ^, ji = — 1, the Hamiltonian 2[C] of reference [13] is 
obtained 

ti 2 -J- t) 2 4- A n 4- fc ~ 16<; - 16 



i? 7 : Class I 3 re = 0, A = 



dinates is given by 



The form of the Hamiltonian in Liouville coor- 



-") V V (-l+ e 2(l+0) 2 (_l +e 2 (»7+€)) 2 



e 2( S -„) e?-" (l+ e 2 (€-<?)) 

(-l+ e 2(?-'7)) 2 (-l+ e 2(«-'7)) 2 

By the coordinate transformation 

1 1 

£ = v + iu, r)=v-iu, p^ = -{p v - ip u ), p v = -{p v + \p u ) 

and putting \i — a, v — 1, the Hamiltonian 4 [5] of reference p3] is 
obtained 

*,»(2«)(rf + rf + lg$ + !S$) +B ,»-m 

# = - 7 r h n 

2 cos(2 w) + a 
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i?§: Class I 3 k = — /i, A = v The form of the Hamiltonian in Liouville co- 



ordinates is given by 

(_ 1+e 2 (>7+0)2 ^i +e 2(l+0) 2 (-l + e2(-7+«)" (_i +e 2 (-„+0 ) 2 



By the coordinate transformation 

£ = arcsinh(tan(0 — i u)), rj = arcsinh(tan(</> + i u)) 
1 



Pt = - cos (<f) - i uj) \Jl + tan 2 (0 - i (p^ + ip w ) 



Pv= 2 cos2 (^ + i u ) V 1 + tan2 (0 + i w ) (P<f> ~ iPu) 

and putting p = a, v = 1, the Hamiltonian 4[C] of reference p3] is 
obtained 

rr P^+pS cos 2 W + cosh 2 M +C3 ^jjggPj^W ) , m+2n-fc-2l 

— " + 2 I Q~ 2 q + 2 1 ct — 2 T 2q 

sinh 2 (2 w) sin 2 (2 <|i) sinh 2 (2 w) sin 2 (2 0) 

Where C = m ( 2 +")+( fc + 2 Q (~ 2 +«)+2 n (2-3 a) £ _ fc(-2+a)+(m+2n) (2+a)-2 1 (2+3 a) 

i = (m+2n) (-2+a)+(fc+2 8 £) (2+a) 8 ° 
3 8 a 

The form of the Hamiltonian in Liouville co- 



R 9 : Class Hi k = 0, A = p 



ordinates is given by 

jj _ PvPs, ^ kr] £ + + mr] + n 

v(v + + » v(v + €) + v 

By the coordinate transformation 

£ = u + if, rj = u - iv, p € = - ip„), p,, = + ip v ) 

and putting // = 1/2, z/ = 0, the Hamiltonian (2) of reference [12] is 
obtained 

Pi + P, 2 , ^ i (£ - m) v k(u 2 + f 2 ) 
# = — — — H h — — + — + £ + m 

4-u u u u 
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Rioi Class Hi A = 0, /i = The form of the Hamiltonian in Liouville co- 



ordinates is given by 

p v p^ n + mr] + + kr]£ 
n = 1 

by the coordinate transformation 

£ = u + iv, rj = u - iv, p 5 = -{p u - ip v ), p v = ^{p u + ip«) 

and putting k = |, v — 1, the Hamiltonian 3 [A] of reference [T3] is 
obtained 



Pu + Pv+^(t + m ) u + 4 i - m) + 4 (n — 4 fc) 



4 + u 2 + v 2 



+ 4k 



Rw'. Class II 2 k — 0, p, — 



ordinates is given by 



if 



The form of the Hamiltonian in Liouville co- 

m+i y/rj {y+Q+k (3 



n + 



+ 



A (77 + + v A (77 + £) + 1/ 

by the coordinate transformation 

f = M + iu, 77 = u - if, p e = hp u - ip v ), p v = hj) u + ip v ) 
and putting A = ^, v = 0, we have 



+ m + (n + 2 £ u) yfu — i v + k (4 -u — 2 i t>) 



4m u a/m — i v 

For this Hamiltonian the additional integrals of motion have the form 



1 k 2 2y/u-\v + nV T r ~ v -ikv] 
A — — X x ^ '- 

2 2 u 



B = Xo 



2v (to (i u+l) +v (zls^II _i fcu)) 



u \/u — iv 
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where K, X%, X 2 are the three integrals of the free motion of reference 



1) / \ 1) / 

K = p v , X 1 = PuPv-— [vl + vl) and X 2 = p v {yp u - up v )-— [p 2 u + p 

This system was not included in reference |12j. Therefore is a new 
superintegrable system, which is studied here for the first time as far 
as we known. 



i?io: Class IIt k 



-A, /j, 



-v 



ordinates is given by 



The form of the Hamiltonian in (£, if) co- 



H 



PiPv+^ + k 4 + ^ + 

+ # + £ + 



n 



By the coordinate transformation 



f = 2 y/uv, r] = i A /- ,P£ 



a pv \ . . f 3/2 . 

Pu + \-Pv) , Prj = "I VUVp u + l—=p v ] 

v V it / \/u 



and putting k = i A = — |,/i = i,i/ = — i the Hamiltonian 3[D] of 
reference [13] is obtained 

m 2 P 2 - ^ 2 P 2 + 2 ( 2 - k ) + i n) u + 2 (2 (£ - k) - i m) v - 2 (fc + 



(u + v)(2 + u-v) 



+2 (Jfe-^) 



i?i3: Class II 3 k = —fi, A = 0, v = The form of the Hamiltonian in (£, 77) 



coordinates is given by 
# = 



P£P„ + + § + S + n 



s 1 i 

f? 3 »7 2 



By the coordinate transformation 



Uy/UV 
U + V 



>uv 
vwuv 



*uv 



1 



P« + - JL —Pv , Pr, = ~7;VUV ( U Pu + VPv) 

u + v 2 



and putting k = —4, /1 = 4, the Hamiltonian 3[C] of reference p3] is 
obtained 



H 



u 2 p 2 u -v 2 p 2 v + ^(u + v) + 2n^-4£^ k 
(u + v)(2 + u-v) 4 
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VIII Superintegrable potentials on a mani- 
fold with curvature zero 



Let us consider the manifold corresponding to the metric in the Liouville 
coordinates: 

ds 2 = g(£,ri) d^dr] 

the curvature is defined by: 

K= ^m iD!i=0 (62) 

The above constraint imposes restrictions on the choice of the parameters 
k, X, \i and v. In Table [3] we can see the possible choices of the values of the 
above parameters. 

In this category there are three families of potentials : 

• The potentials on the complex E plane corresponding to the Hamilto- 
nian: 

H = P 2 X +P 2 y + V(x,y) 

these potentials are classified in reference [5] and finally the final list 
of potentials are given in reference [6]. In this section we follow the 
enumeration of the potentials given by this exhaustive list [6] 

• The Drach potentials corresponding to the Hamiltonian 

H = VxVy + V(x, y) 

these potentials are classified in reference [8]. In this section we fol- 
low the enumeration of this reference for the Drach potentials with 
quadratic integrals of motion. 

• The potentials on the real hyperbolic plane H2 corresponding to the 
Hamiltonian: 

H = p 2 x -p 2 y + V(x,y) 
these potentials are classified in reference [9] 

Generally the condition fl62l) restricts the choices of the constants k, A, \i and 
v in one independent parameter. These values characterize the permitted 
metrics in the classification given in Section |VJ The permitted choices are 
given by the following list 
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Fx: Class I x k = 0, A = 0, fi = 



coordinates is given by 



The form of the Hamiltonian in Liouville 



„ p v P ( k(Z + v) + 4l(Z + v) 2 -nt-v) 2 + T^+n 

H = 1 

v v 

By the coordinate transformation: 

£ = x + i y, rj = x - i y, p^ = ^ (p x - ip y ), p v = ^ (p x + ip y ) 

and putting v = ~, the complex plane Hamiltonian Ei of reference [6] 
is obtained 

H = P 2 X +P 2 y + 16 £ (4x 2 + y 2 ) + 8 kx - — +4n 

Also the potential V a of reference [9]is generated. 

By the coordinate transformation: 

. x . 
q = —, i] = y, p^ = rp x , p v = p y where r is a constant 

and putting v = r, the Drach Hamiltonian (a) (r 7^ 0) of reference [8] 
is obtained: 



k(x + ry) mr 3£(x 2 + r 2 y 2 ) 10£xy n 

r 2 (x — ry) 2 r 3 r 2 r 



This system is generated by the two systems of revolution Ri and R2 
see Table |6j In these cases two parameters among the k, A, ft and v are 
zero. When a third parameter is zero these systems are degenerated to 
the system F\ in Table [3j 



F 2 : Class I 2 k = 0, A = 0, fi = 



The form of the Hamiltonian in Liou- 



ville coordinates is given by 



v v 
By the coordinate transformation: 

i = x + i y, if = x - i y, p^ = - (p x - ip y ), Pr, = \ (Px + [ Py) 
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and putting v = ~, the complex plane Hamiltonian E\ of reference j6] 
and the Hamiltonian V b of reference [9] are obtained 

H = P 2 X + P 2 y + 16 £ (x 2 + y 2 ) + 4 - + 4n 



By the coordinate transformation: 



£ = — , rj = y, = rp x , p v = p y where r is a constant 

and putting v = r, the Drach Hamiltonian (b)(r ^ 0) of reference [8] is 
obtained: 



Aixy mr kr n 
H= PxPy + — ^- + -2 + — — -2 + - 

t [x — ry) (x + ry) r 



F 3 : Class I 2 k = 0,;U = 0, z/ = The form of the Hamiltonian in Liouville 



coordinates is given by 

4A£r? 4A£r/ 
By the coordinate transformation: 



x + iy x-iy 

Pi = y/2(x + iy) (p x - \p y ) , p v = ^2(x -ly) (p x + \p y ) 
and putting A = 1, the complex plane Hamiltonian £7 16 of reference j6] 

H = pl + p 2 y + - t 4=^2 f? + I + nm ) +i 

y \*x z + y z \2 x + v x z + y z x — +y / 
and the Hamiltonian V c of reference [9] 

n k+m k-~m 



£[ _ p2 _|_ ^2 _| 2 | 2 _| 2_ 

x y a/x 2 + y 2 y 2 ^fxP^Ty 2 y 2 



are obtained 
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By the coordinate transformation: 

£ = y-,v = Vy,Pt = zVrxpx , Pr? = 2y/yp y 

and putting A = r, the Drach Hamiltonian (g)(r ^ 0) of reference [8] 
is obtained: 



# = PxPy + —p= + 



(m—k) 



k+m 



+ 



F 4 : Class Hi A = 0, li = 0, v = 



a^y (x — ry) yjxy (x — ry) r 

The form of the Hamiltonian in Liou- 



ville coordinates is given by 

p v p$ n + mr] + (£ + krj) £ 
ti = H 

By the coordinate transformation: 



\J x+yfx 2 +y 2 +yj x- y/x^+y 2 yj x+y/~x 2 +y 2 —yj x-y/xP^y 

* = 73 ' ^ = 73 



\/ x +V x 2 +y 2 + \ x-J x 2 +y 2 

Pi = — 7T^ P-+ 



: T^/ x ~\/ x2 +y 2 + \/ x +\f e 2 +J/ 2 ^ +\/ x 2 +J/ 2 x~y/x 2 +y 2 + ^J x+y/x 2 +y' 



V2y/x 2 +y 2 



V2i 



and putting k = 1, the complex plane Hamiltonian i?20 of reference j6] 
and the Hamiltonian V d of reference [9] are obtained 



:Py 



H = p 2 x +P 2 y +7=(n++^V x + V?+? + 



a/x 2 + y' 2 ) + A; 



x ^ yj x 2 +y 2 

By the coordinate transformation: 

£ = 2y/x,r) = 2y/y,p( = \fxp x , p v = y/yp y 

and putting k = |, the Drach Hamiltonian (c) of reference [8| is ob- 
tained: 

n 2m 2£ 
H = p x p-y -\ — — = H — -j= H — — + Ak 
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F 5 : Class 1^ « = 0, A = 0, v = 



ville coordinates is given by 



The form of the Hamiltonian in Liou- 



jt _ p v Pz + krj£ + + mr] + n 



By the coordinate transformation: 



£ = x + i y = z, rj = 2 Jx — i y = 2\fz 



Pi = ^(p x - ^ Py) , Pr, = ^yx - i y (p x + i p y ) 

and putting fj, — |, the complex plane Hamiltonian i? n of reference jH] 
is obtained 



2! 4n 

Using the coordinate transformation: 

£ = Xj r) = 2 y/y, p£ = p x , p v = jyp y 

and putting fi = ~, the Drach Hamiltonian (e)(r = 0) of reference [8] 
is obtained: 

# = PzPu H ^2kx + £ h 2m 



F 6 : Class II 2 k = 0, A = 0, ^ = 
ville coordinates is given by 



The form of the Hamiltonian in Liou- 



Pn Ve k£ ££ 3 k ,/rj £ri m n 
v v v v \p] v v 

By the coordinate transformation: 

£ = x + iy, ri = x-iy,p£ = ^{p x - i p y ) , p v = ^{p x + i Py) 

and putting u = j, the complex plane Hamiltonian £"9 of reference [6] 
is obtained 



uk (2x — i y) . 4m 

H = P 2 x + pI + , ^ +8^x + ^= + 4n 

y/x-iy \/x - 1 y 
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By the coordinate transformation: 



x 



£ = - , V = V , = r Px , Pr, = Py 



r 



and putting u — r, the Drach Hamiltonian (e)(r 7^ 0) of reference [8] is 
obtained: 



I (x + 3 ry) t n 



H = p x p y + + — (x + ry) + — h 



F 7 : Class II 2 k = 0, A = 0, v = 



ville coordinates is given by 



The form of the Hamiltonian in Liou- 



rr PvPtVV , k £ , ^V^^ , 3k V , ^fa)* , m , 

-n = 1 1 1 1 h 1 

jj n n fi fi fi fi 

Using the coordinate transformation: 

(x — i y) 2 ~z 2 



£ = x + i y = z, 77 



2 



Pi = \ (Px ~ i Py) , Pv = ~ 2 (x-iy) ^ ~ 1 P ^ 

and putting /x = , the complex plane Hamiltonian E 1Q of reference 

[6] is obtained 

3 1 

H = P 2 X +P 2 y - 4nz + 4 v^i fe (2 - - z 2 ) - 4£(zz - - z 3 ) + 4i V2 m 

By the coordinate transformation: 

£ = x , 77 = ?/ 2 , p € = p x , ^ = ^p y 

and putting /j = ~, the corresponding Drach Hamiltonian is obtained: 

H =p x p y + 2k (x + 3 y 2 ) + 2 1 y (x + y 2 ) + 2 n y + 2 ■ 
This potential was not included in the list of reference [S|. 
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F 8 : Class II 3 k = 0, A = 



ordinates is given by 



The form of the Hamiltonian in Liouville co- 



rr PvPi , k £ , . m n 

ti = 1 — 1 1 — - 1 — 

v rf v v r\ A vv 

Using the coordinate transformation: 

£ = x + i y = z , t) = x - i y = z , p^ = - (p x - i p y ) , p v = - (p x + i p y ) 

and putting /i = , v = | , the complex plane Hamiltonian E 8 of 
reference [6] is obtained 

2 2 4m 

# = Px + Py + ~=3~ + ~=2- + ^^ZZ + An 

By the coordinate transformation: 

^ = x, r] = y,p^ = p x ,p v =p y 

and putting v — 1, the Drach Hamiltonians (/) of reference [8] is 
obtained: 

m kx 
H = p x p y + Zxy+^^ + n 



F q : Class II, k = 0, A = 



nates is given by 



The form of the Hamiltonian in (£, 77) coordi- 



P,P, , + 



By the coordinate transformation 



1 _ x - i y - J(x - iy) 2 + 16 //z/ 
£=-(x + iy),7i = 



2 v (p x + i p y ) J(x - i |/) 2 + 16 (i v 

Pi = Px ~ 1 Py , P v = - 



x — iy — J(x — i y) 2 + 16 // z/ 
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and putting // = — |, v — | the complex plane Hamiltonian £7 of 
reference [6] is obtained 



\J w 2 -c 2 ^w 2 -c 2 (w+y/w 2 ^ 2 ^ 



4fc — 1 2(m—n) 
— W W ~\ * '- 

c 

where w = x + iy,W = x — iy 
Using the coordinate transformation: 

X + y/x 2 + 4/iZ/ 2 z/ v / x 2 ~+T/7z/ 

and putting z/ = ^ , /1 = ^ the Drach Hamiltonian (i) (r 7^ 0) of 
reference [8] is obtained 

H = p x p y + xy + 7== + r^— — 

The form of the Hamiltonian in Liouville coor- 



F l0 Class II 3 fi = 0, 1/ = 

dinates is given by 



Using the coordinate transformation: 



£ = 2Jx-i y = 2v^, 77 = 2 a/x + i 2/ = 2x/i 



and putting k — , A = ^ , the complex plane Hamiltonian £^7 of 
reference [5] is obtained 

9 9 k „ m An 

H = P x + P y + — + 16 1 + —j= + 



2 zyzz y zz 

Using the coordinate transformation: 

£ = 2 Vv, V = 2 y/x, Ps = y/ypy , Pr, = \fxp a 
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and putting k = 0, A = 7 , the Drach Hamiltonians (d) of reference jH] 
is obtained: 

k „ m n 
H = p x Py + - — g + 4 £ + r + 



4x 2 4^12 v^v 7 ^ 



nates is given by 



Fn: Class II3 fi = 0, z/ = The form of the Hamiltonian in (£, 77) coordi- 



7J pep, h (»+***) + 



?? 



By the coordinate transformation 



£ = 2 Jx + i y, 7] = \ fx - i y + J (x - i y) 2 - 4 



1 , . v / — (Px + i Pi/) A/(x-i?/) 2 -4 

P( = 2 (P* ~ 1 P» W x + 1 y ' ^ = = 

yx -i y + ^(x - i y) 2 - 4 

and putting k — — | , A = | we have the Hamiltonian in the complex 
plane E%g of reference [B] 

By the coordinate transformation: 



\/x 2 — r 2 



i = 2^/y ,ri = 2\lx + Vx 2 - r 2 , p$ = ^/y p y , p v = = Pr 

+ \J x 2 — r 2 



and putting k = — 2r 2 , A = | the Drach Hamiltonian (/i)(r 7^ 0) of 
reference [8] is obtained 

4nr + m 4nr — m ( 16 i r 2 +k\ ^ n 9 . 

^ 472 J # lb-Cr^ — ft 



# - P*P y + , + I + —7=1==^ + 

Jy(x — r) Jy(x + r) vx —r 
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IX Superintegrable potentials on a manifold 
with constant curvature 



Let us consider the manifold corresponding to the metric in the Liouville 
coordinates: 

ds 2 = g(£,ri) d^dr] 

the curvature is defined by: 

1 d 2 

K = — — — — In g = Constant (63) 

2g o^oi] 

The above constraint imposes restrictions on the choice of the parameters 
k, X, n and v. In Tabled we can see the possible choices of the values of the 
above parameters. 



d: Class Ii k = 0, A = 0, fi = 1/K, u = 



The form of the Hamiltonian 



in Liouville coordinates is given by 



For K — 1 and using the coordinate transformation: 

cot(|) (— i p^+pe sine) 



£ = e^tan(§), p e ^ 

r, = -j* cot(f ) , ^ = ^!!%H 

the Hamiltonian S\ of reference [6] in spherical coordinates 9, <fi is ob- 
tained: 

•> Pi , <k,*cos0 „ A -, (1 + 2 cos(2fl)) e 2i(t> 
sm v sm V sm V sm v 



C 2 : Class I 2 k = -1/K, A = 0, /i = 0, z/ = 



The form of the Hamilto- 



nian in Liouville coordinates is given by 
H= -Kp v pz(r] + £) 2 - 
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Class 


K 


A 




V 


Plane 
potentials 
Irom ret [oj 


Drach potentials 
from ref [8] 


Potentials 
from ref [9] 




h 













E 2 


(a) (r ^ 0) 




F 2 


h 















(6) (r ^ 0) 




F 3 















Exe 


(g) (r + 0) 


V c 


F 4 


Ih 













E20 


(c) 


yd 


F 5 















En 


(e) (r = 0) 




F G 


Ih 













E 9 


(e) (r 0) 




F 7 















Exo 






F s 


Ih 












E s 


(f) 




F 9 














E 7 


(i) (r ^ 0) 




Fio 














En 


(d) 




F u 














Exg 


00 (r 7^ 0) 





Table 3: Potentials with curvature zero and two quadratic integrals of motion 





Class 


K 


A 


n 


V 


Potentials 
on the sphere (K=l) 
from Ref [6] 


Potentials 
from Ref [9] 


Cx 


h 








1/K 





Si 




c 2 


h 


-1/K 











s 2 




c 3 




-1/K 





1/K 





S4 




c 4 


h 








A/K 







U c 


c 5 










2/K 


1/K 


Sg 


u a ,u b 


c Q 




-2/K 


-1/K 


2/K 


-1/K 


s 7 


u e 


C 7 




-A/K 





A/K 





S 8 





Table 4: Potentials with constant curvature and two quadratic integrals of mo- 
tion 
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For K = 1 and using the coordinate transformation: 



£ = -±ie^cot(§), p c 



tan(|) (p<j,— i pe sine) 



^-lie^tand), ^ = 

the Hamiltonian of reference [6j in spherical coordinates 9, <fi is ob- 
tained: 

H = pj + - A;-£e 4i0 esc 2 9 - m sec 2 + n e 2i esc 2 9 

sin 6* 



C 3 : Class I 2 « = -1/K, A = 0, /i = 1/K, v = 



nian in Liouville coordinates is given by 



The form of the Hamilto- 



H K . P., PL_ , K + £ + «! - f « - -rt 2 + ggg + n 
For .fT = 1 and using the coordinate transformation: 



£ = 2 J-i e» * cot(f) , -;^ m " 

V V2/ ' -fS 2 ->/ — i e 1 <t> cot(|) 

V -Zy/ie tan( 2 j, p, - 2 ^ - - 

the Hamiltonian S4 of reference [5] in spherical coordinates 9, <fi is ob- 
tained: 

H = pl + + (m-k) - i (jfc + m) cot# - 64£^4 - 4n-^- 

shr6> sm z 9 sm. 9 



C 4 : Class I 3 k = 0, A = 0, fi = A/K, u = 



Liouville coordinates is given by 



M — 4e 2 (-i+€) + 



The form of the Hamiltonian in 



+ 



(-l+e 2 (-v+0) 2 K [ '— + L 



4 e 2(-7j+0 

For .fT = 1 and using the coordinate transformation 
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£ = \ ln(-i e 1 cot(|)) , p ? = -i ^ - p e sin 9 
7] = | ln(i e 1 * tan(|)) , p,, = -i + p<9 sin 

the above Hamiltonian in spherical coordinates 9, <fi has the form 

TT o Pd> k £cos(j) m ni 

H = Vn H + o— g 1 5— q 1 cot 6* 

6 sin 2 9 4 sin 2 sin 2 2 sin 2 9 sin 2 4 2 

This is the Hamiltonian U c of Ref. [S] in spherical coordinates. This 
potential is missing in the classification scheme of Ref. [6]. 



C 5 : Class I 3 « = 0, A = 0, /i = 2/K,v= l/K 



nian in Liouville coordinates is given by 



The form of the Hamilto- 



H K 2e *«-,) P + " e«-,( 1+e2 + 

(_l + e 2 (€-»7)) 2 (-l + e 2(C-'))) 2 

e 2(, + C) fe +e "+«(l +e 2 ^))^ + e 2(^) m e-»»+t(l +e »tt-^)n 

(_i +e 2 („+o) 2 (_i +e 2 (v+^y (_i +e 2«-^)) 2 (_i +e 2 ci^y 2 

^2^0 + ^-"(l+« 2 «"" ) ) 
(_l +e 2 C«-^)) 2 (-l + e 2(«-»?)) 2 

For K = 1 and using the coordinate transformation: 

£ = ln(e^tan(§)), Pi = sing 
77 = ln(-e^cot(|)), P<7 = - i^±^ 

the Hamiltonian 5g of reference j6] and the Hamiltonian £/ b of ref. j9] 
in spherical coordinates 9, <fi are obtained: 

H = p] + + 

+ + (jfc + 2 £ cos(2 v?)) esc 2 9 csc 2 (2 0) - sec 2 

while using the coordinate transformation: 

t -i ,1 /l+i e 1 ^ tan(f)x , , 1-i e 1 * cot(| ) x 

e= _^ + l n (___|_) ) ^ = _ ff + ln( __2 )? 
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-(p (cot(f)+e 2i *> tan(f))) |pecos(|) 2 (l+e 2 1 ^ tan(§) 2 ) 
V(, 4 e i f ~^~4> i 

^Ve (l+e 2i ^+(-l+e 2i <t>) cos(0)) p0 (e 2i * cot(|)+tan(|)) 

we have the Hamiltonian f/ a of ref. [9] in spherical coordinates. 
H = Pe + 3e+ 2 2 

I m+2n _|_ (2n-m) esc 2 6* esc 2 , k-2l , k+2 I 

4 4 2(cos6»-cos0 sing) 2 2 (cos 6»+cos </> sin6») 2 



C 6 : Class I 3 k = -2/K, X = -1/K, \i = 2/K, u = -l/K 



The form of the 



Hamiltonian in Liouville coordinates is given by 



H= K 



(_ 1 + e 2 (S-V)y (_l + e 2(£-„)) 2 (-l + ^+O)* (_ 1 + e 2 („ + «)) 2 

e2(r; + e)fe +e ^(l + C 2 ^ + 0) ;+ eH - n+Qm + e«-'<(l + e 2 «-!))„ 
(-l+eZ^+C)) 2 (_l +e 2(^+«)) 2 ^i +e 2(«--,)) 2 (_ 1+e 2(e-^)) 2 
+ 2e 2( 5 -,) _ e ^(l + e 2 («-' ) ))_ 2e2( „ + 

(_l +e 2 K^)) 2 (, 1+e 2( € -^)) 2 (_ 1+e 2 („+?)) 2 (_ 1+e 2(r,-R)) 2 

For .fT = 1 and using the coordinate transformation: 

^axcsinh(-ie>«cot(§)), P{ = V^^~5> 

77 = arcsinh(-i e 1 * tan(f )) , p v = ^ - 2 1 t L(|) ^ — ~ 

the Hamiltonian SV of reference j6] and the Hamiltonian U e of ref. [9] 
in spherical coordinates 9, <p are obtained: 

£7 _ „2 1 p l I -fc-2l+m-2n , (fc-2 l-m-2 n) sec(6>) 2 . 

^ - + s + 8 + 

I (fc+2Z+m-2n) cos(<ft) sin(6>) (fc-2 Z+m+2 n) scc(6>) sin(</>) tan(fl) 
v /cos(6») 2 +sin(6») 2 sin(0) 2 8 v /cos(6») 2 +sin(6») 2 sin(</>) 2 

The form of the Hamil- 



C 7 : Class I 3 « = -4/iT, A = 0, // = i/K, u = 



tonian in Liouville coordinates is given by 

(_l+ e 2 («-.?)) 2 (_i +e 2 (>7+C)) 2 

1 e 2 fajj) > , + C ( 1 + e2 ( " + ) ' , e 2 «Z2) m , ^" ( 1 + e2 ^1) 

+ K 



e 2 («-»?) 4 c 2 (1+0 



(_l +e 2(C-))) 2 (_i +e 2(-7+0)' 
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For K = 1 and using the coordinate transformation 



1 (i (-l + Jl + (a + ir) 2 )\ 1 /i (1 + + r ) 2 



£ = o ln n ' *7 = o ln 



2 \ j + i r / 2 \ cr — it 



^l + ^ + ir) 2 ((o- + ir)p CT + (-i o- + r)p T ) 



p,, = i ^/l + (cr - i r) 2 ((i a + r) p CT - (cr - i r) p T ) 

the Hamiltonian Ss of reference [6] in horospherical coordinates cr, r is 
obtained. 

, 2 2 \ 2 (fc+m) (-1+ct 2 +t 2 ) 

V 7 8 v^ 4 + (-1+t 2 ) 2 +2<7 2 (1+t 2 ) 
+ (n = ^£ = r+r 2 ) (l+n) ( CT 2 +r+r 2 ) + m _ fc 

4y / (<7 2 +r 2 )( ( 7 2 + (-l+T) 2 ) 4 y / ( CT 2 +r 2 ) ( ( j 2 + (1+t) 2 ) 8 

We use the inverse transformation 

2„ j „2e _ R 2„ -2{ 



i _ 1 e 



l+e 41 ? l+e 4 «' l+e 4 'i l+e 4 « 

for the verification of Poisson brackets. The above transformations are more 
appropriate for the corresponding calculations 



X Superintegrable systems with a linear and 
a quadratic integral 

In the case of a linear integral of motion and a quadratic integral of motion, 
there is a Liouville coordinate system where the Hamiltonian and the linear 
integral of motion are written as: 

I = (P« + Pvf or A = (P* - Pvf 

and 

„ = PtPv , 9(£-V) R= PjPy , f(£ + V) 

c(e - v) ^(e - v) + v) + v) 

From the forms of the integral of motion, which are given in Section |V] we 
can find all the possible subclasses corresponding to a linear and quadratic 
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integral of motion in Liouville coordinates. In all the above cases we remark 
that the potential depends on two parameters among the k, £, m and n. In 
Table [5] we give the possible cases of superintegrable systems with a linear 
and a quadratic integral of motion: 

In Table [6], the possible superintegrable systems which are defined on a 
surface of revolution with a linear and a quadratic integral of motion are 
listed. 

In Table [7J the possible superintegrable systems which are defined on a 
surface of zero curvature with a linear and a quadratic integral of motion are 
listed. 

In Table [HI the possible superintegrable systems which are defined on a 
surface of constant curvature with a linear and a quadratic integral of motion 
are listed. 

XI Discussion 

The findings of this paper are summarized as follows: 

1. The superintegrable systems with quadratic integrals of motion can be 
classified in six subclasses. Each subclass depends on seven parameters. 
Four among these parameters («, A, fi and v) determine the metric of 
the manifold, on which the system is determined. These parameters 
are characteristic of the system's "kinetic" energy. The remaining three 
parameters define the potential (The potential depend on four param- 
eters k, £, m and n but only three among them are independent). 
For each subclass, the analytic explicit forms of the manifold metric, 
the potential and the integrals of motion are calculated. Also the con- 
stants of the quadratic Poisson algebra of integrals of motion are given 
as functions of the energy and the eight parameters are given. 

2. All the known two dimensional superintegrable systems are systems de- 
fined on surfaces of constant curvature or on surfaces of revolution. All 
these systems are classified in these six classes. Each class is character- 
ized by the values four parameters k, A, /x and u, which are determined 
by the form of the assumed manifold. If we fix the manifold, let us 
suppose that the manifold is a manifold with negative constant cur- 
vature, the possible values of the of the parameters k, A, fi and v are 
calculated for each subclass. Therefore we can "guess" the existence 
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Class 


K 


A 


fi 


V 


k 


£ 


m 


n 


GLi 


h 






















GL 2 
























GLx) 


h 






















GL 3 
























GL4 
























GL 5 


h 






















GLi) 




k = ^2A 




H = ±2u 


k = t2£ 







m = ±2n 




GLq 


Ih 





A = ±/i 









I = ±ra 






GL 7 


Ih 






















GL 3 ) 


Ih 























Table 5: General superintegrable integrable systems with a linear and a quadratic 
integral of motion 























Potentials 


Potentials 




Class 


K 


A 




V 


k 


/ 


m 


n 


by revolution 


by revolution 






















from ref [TBI 


from ref [12] 


RLi 


h 






















2.2[D] 




RL 2 


























(3) 


RL 3 


h 






















A.2[D] 




RL^ 


Ih 






















3.2[E] 





Table 6: Potentials by revolution with a linear and a quadratic integral of motion 
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Class 


K 


A 


n 


V 


k 




m 


n 


potentials 
from ref [6] 


FLi 


h 























E 6 


FL 2 

























E 5 


FL 3 


h 























E 3 


FLi 



























FL 5 


Hi 























E 4 


FLq 


lh 























E13 


FL 7 


Ih 























E u 


FL 8 
























E\2 



Table 7: Potentials with curvature zero and with a linear and a quadratic integral 
of motion 



Class 


K 


A 


n 


V 


k 


/ 


m 


n 


Sphere 
potentials 

(K=l) 
from ref [6] 


h 








l/K 

















h 


-l/K 





















5 3 




-l/K 





l/K 















Sq 



Table 8: Potentials with constant curvature and with a linear and a quadratic 
integral of motion 
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of the permitted superintegrable systems and to classify these systems 
in tables. Using this technique we can classify all the possible known 
superintegrable systems and to investigate the possible missing poten- 
tials. With this method a new superintegrable system was found for 
the class I superintegrable systems by revolution, given in ref [5] . 

3. Generally for any values of the parameters k, A, \x and v the associ- 
ated manifolds are neither surfaces of constant curvature nor surfaces 
of revolution. Therefore we have investigated superintegrable systems, 
which are not yet known. We believe that all the possible superinte- 
grable systems with two quadratic integrals of motion are investigated. 

4. The two dimensional 'non degenerate' [5] superintegrable systems are 
classified by the values of the constants of the Darboux equations [9] [5] 
and the constants of the system. The relation of these constants with 
the constants of the quadratic Poisson algebra is explained. 

5. The six classes of superintegrable systems are the equivalence classes 
of Stackel equivalent systems. 

From the above discussion several open problems arise: 

• The superintegrable systems for the case of cubic integrals of motion 
are under investigation [SJ [31], [32] [33] [M] [35] [36] . The general structure 
of these systems is recently investigated [H] but the general form and 
their classification is not yet known for manifolds which carry integrable 
systems with one third order integral of motion. 

• The quantum counterparts of the general six subclasses of superinte- 
grable systems with quadratic integrals of motion are not yet known. 
In Section [III the separation of variables of these systems has been 
explicitly written. The form of the separation of Schroedinger equa- 
tion can be written in a Liouville coordinate system. The solutions 
of the quantum Schroedinger equation can be calculated. This work 
is under current investigation. The form of the Poisson algebra can 
be generalized in a quadratic associative algebra, whose energy eigen- 
values are generally calculated by using deformed oscillator techniques 
[T8] [T9| 1201 I2TI [22] . From the form of the Poisson algebra, one can be 
guess that the energy eigenvalues of these quantum systems are roots 
of cubic polynomials. 
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• The general form of the three dimensional superintegrable systems with 
quadratic integrals of motion is not yet known. The Poisson algebras 
and the associated quantum counterparts for the three dimensional 
superintegrable systems are not yet fully studied. Recently [37J EH] 
the quantum quadratic algebras have been written down, which are 
not generally closed as polynomial algebras. A systematic calculation 
of energy eigenvalues with algebraic methods has not been performed 
yet. 
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Appendix 



A Polynomial combinations of integrals 

Let consider a superintegrable system with two quadratic integrals of motion. 
The general forms of the Hamiltonian and the integrals of motion in Liouville 
coordinates are: 



A = p\ + kp%- 2c(£, V)j§r% + , V), k = or 1 (64) 

5 = « 2 (0p| + b 2 (v)p 2 v - Mt,v)ffi% + Q(Z,V) 

In this Appendix we consider that the system has quadratic integrals of 
motion. We assume that the systems has not any linear integral of motion. 
That assumption excludes the system H = p\ + Py, because it possesses two 
linear integrals of motion. In this Appendix, we shall prove the following 
proposition: 

Proposition 4 Let M be an integral of motion, which is a polynomial func- 
tion of the momenta of even order. We assume that this integral contains 
only monomials of momenta of even order, i. e. 

2n 

M= £ a M (e,i7)pJp{ (65) 
k+e=even 

Then M is a polynomial of order n of the integrals H, A, B . 

The existence of three integrals of motion implies that the integral M is 
functionally dependent integral, i.e. there is some smooth function $ (non 
generally a polynomial one) such that: 

$(M, A, B, H) = or M = f(A,B,H) 
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but the function f(A, B, H) is not evident that it is a polynomial one. From 

A + 2cH -q-2cV k 
B + 2(3H -Q-2/3V b 2 (^7]) 
1 ~k I 

1 A + 2cH -q-2cV 
a 2 fay) B + 2{3H - Q - 2(5V 

I * * I 

H-V 

The above equations imply that the momentum monomials p| p e v inside the 
sum sign in (1651) can be written as polynomials of the integrals A,B,H with 
coefficients, which depend on £ and rj. Therefore the integral of motion M is 
written: 

M= £ c ijk (^rj)A i B^H k 

0<i+j+k<n 

Generally the coefficients should be constants not depending on the vari- 
ables £,7/. If these functions are non constant functions %■&(£, rj) then we 
choose a fixed value of these parameters £o, T]q. In general there is an infinity 
of trajectories passing through these position values £o,?7o- Each trajectory 
is characterized by a special value of the integrals H, A, B, therefore M is 
a polynomial of fixed constants for an infinity of trajectories. For another 
pair of parameters £1,771 we have another choice of the coefficients in (|65l) 
therefore there is a relation of the form: 

E m) - Cijkfa, r/o)) A i B^E h = 

0<i+j+k<n 

The above relations means that the integrals of motion H, A, B are not func- 
tionally independent functions, that is a contradiction to assumption initial 
regarding the independence of the integrals. Therefore we have proved that 
the polynomial expansion (1651) of the integral M is indeed unique, when M 
is a even polynomial of the momenta. 

A direct application of Proposition H] is the following Proposition: 

Proposition 5 Let a superintegrable two dimensional system have even qua- 
dratic integrals of motion H,A,B. If we put C = {A,B}. The integrals 



we can see that 



PtPr, = g 
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{A, C} and {B,C} can written as quadratic polynomials of the integrals. 
The integral C 2 is a cubic polynomial of the integrals H,A,B. 

The above proposition was taken as a conjecture in the previous work 
[32] [22] an d [23] [2S]- Here we prove that this assumption is indeed true. 
A generalization of Proposition [5] is indeed true for superintegrable two di- 
mensional systems with an integral, which is an odd cubic polynomial in 
momenta [32| |3T] . This Proposition means that the superintegrable two di- 
mensional systems with even quadratic integrals correspond to a quadratic 
Poisson algebra, which is characteristic for the superintegrable system. The 
situation in three dimensional superintegrable systems is not clear J37J, EB] 

Acknowledgements: One of the authors (C. D.) would thank W. Miller 
and J.Kress for their valuable comments and for indicating us the substantial 
references [14,15,16], which treat the problem of the classification from the 
point of view of equivalence classes of Stackel transform. 
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